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Abstract 

We construct two anomaly-free supersymmetric U(l)' models with a se- 
cluded C/(l) / -breaking sector. For the one with Eq embedding we show that 
£T) \ there exists a strong enough first order electroweak phase transition for elec- 

OS ■ troweak baryogenesis (EWBG) because of the large soft trilinear terms in the 

Higgs potential. Unlike the Minimal Supersymmetric Standard Model (MSSM), 
the stop masses can be very heavy. We then discuss possible CP violation in 

T 1 

| the Higgs sector, which can be both spontaneous and explicit, even at tree 

level. The spontaneous violation provides a direct source for baryogenesis, 
while its magnitude is mediated by an explicit phase from the secluded sector. 
These new CP sources do not introduce significant new contributions to electric 
dipole moments. EWBG in the thin wall (r leptons) and thick wall regimes 
(top squarks, charginos and top quarks) are systematically discussed. We find 
that the CP- violating stop and chargino currents are very different from those 
obtained in the MSSM. Due to the space-dependence of the relevant CP phases, 
they do not require a variation of tan (3 in the bubble wall to have a non-trivial 
structure at the lowest order of Higgs insertion. In addition to r leptons, top 
squarks and charginos, we find that top quarks can also play a significant role. 
Numerical results show that the baryon asymmetry is large enough to explain 
the cosmological observation today. We illustrate that EWBG and neutralino 
cold dark matter can be accommodated in the same framework, i.e., there ex- 
ists parameter space where a strong enough first order EWPT, large CP phase 
variations across the bubble wall, a reasonable baryon asymmetry as well as an 
acceptable neutralino dark matter relic density can be achieved simultaneously. 
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1 Introduction 



The baryon asymmetry of the Universe (BAU) has been precisely measured by WMAP. 
Combining their five year results with other CMB and large scale structure results, 
they obtain [TJ 

tt b h 2 = 0.02265 ± 0.00059 , (1) 

in excellent agreement with the 95% CL. range 0.017 — 0.024 obtained from big bang 
nucleosynthesis [2]. This implies 

(6.21 ±0.16) x 10" 10 , (2) 
(8.82 ± 0.23) x 10~ n , (3) 

for the ratios of baryon density (%) to photon density (n 7 ) and entropy (s). 

To generate the baryon asymmetry, the Sakharov criteria [3] must be satisfied: 
(1) Baryon number (B) violation; (2) C and CP violation; (3) a departure from 
thermal equilibrium (or CPT violation). The first criterion is obvious. The second 
is required because if C and CP are exact symmetries, no net baryon number can be 
produced. Moreover, the equilibrium average of B is zero in thermal equilibrium by 
CPT, so the third criterion is necessary. 

There are several baryogenesis scenarios: Grand Unified Theory (GUT) baryo- 
genesis [U EH |6] , Affleck-Dine baryogenesis [7] , leptogenesis [H [9] , electroweak baryoge- 
nesis (EWBG) PZH Ellin], etc. (For reviews, see [131 El US US]-) We will concentrate 
on EWBG, where the Sakharov criteria can be satisfied. Baryon number violation 
is induced by the anomalous sphaleron process, which violates B + L but preserves 
B — L at temperatures higher than the critical temperature T c of the electroweak 
(EW) phase transition (EWPT) [TTJ, [181 HHj • I n addition, C is violated because the 
left- and right-handed fermions transform differently under SU(2)l x U(1)y- CP 
can be violated in the Standard Model (SM) and Minimal Supersymmetric Standard 
Model (MSSM) and their extensions by the Yukawa interactions, \x term, and soft 
supersymmetry breaking terms. Non-equilibrium processes may be obtained through 
the expansion of bubbles containing the true vacuum after the EWPT. 

However, to preserve the generated baryon asymmetry after the EWPT, one 
needs a strong first order EWPT so that the sphaleron effects in the bubbles of true 
vacuum are shut down quickly. This requirement can be expressed as3 

^>l, (4) 

c 

where in the SM v(T) is the vacuum expectation value of the Higgs field H° at 
temperature T and T c is the critical temperature of the EWPT. Similarly, in the 

1 Some authors impose the weaker constraint y/2v(T c )/T c > 1. 
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MSSM 



v(T) = a/ t>i(T) 2 + v 2 {T) 2 , (5) 

where v±(T) and VziT) are respectively the VEVs of the two neutral Higgs fields H® 
and H®. Our normalization is such that v (0) ~ 174 GeV. 

In the SM, the EWPT can be strongly first order only if the Higgs mass is 
lighter than 40 GeV, which is inconsistent with the lower bound of 114.4 GeV from 
LEP [20] (recently, it was noticed that this tension can be relaxed in the SM, with an 
extra singlet introduced [2T]). Furthermore, the CP violation from the CKM matrix is 
too small to generate a sufficient baryon asymmetry during the EWPT. In the MSSM 
there may be a small parameter space for EWBG [22J-J4T]. There are additional 
sources of CP violation associated with the phases of the \i and supersymmetry 
breaking parameters. However, a strong enough first order EWPT requires that the 
lightest stop quark mass be smaller than the top quark mass ~ 173 GeV, called 
the light stop scenario [2H [291 EQj ED E2] (for recent work, also see [331 El [36j 
1371 1381 [39]). Also, the mass of the lightest CP-even Higgs must be smaller than 
125 GeV, which leaves a small window above the current limit. In the Next-to- 
Minimal Supersymmetric Standard Model (NMSSM) [42] there is a trilinear term 
A h hSH u H d in the Higgs potential, where S is a singlet under the SM gauge group. 
This allows a strong enough first order EWPT [431 I4"41 14"5*| I4"6] . In the best-motivated 
versions of the NMSSM (motivated by the fi problem [47]), the effective /i parameter 
is given by h(S) (for recent studies, see [H]). However, the original versions involve a 
discrete Z3 symmetry and serious cosmological domain wall problems [49] . These can 
be avoided by introducing an elementary bilinear /i term in the superpotential (or 
linear or quadratic terms in S), but at the cost of reintroducing the /i (or analogous) 
problems [441 143] . The domain wall problem can also be avoided in an alternative 
version of the NMSSM [5U], and in the nearly Minimal Supersymmetric Standard 
Model (nMSSM) [STJ [52J [53] , using a different discrete symmetry broken by a loop- 
induced tadpole. (For reviews, see [54].) 

On the other hand, the possibility of an extra U(l)' gauge symmetry [55] 
is well- motivated in superstring constructions [56], in grand unified theories [57J, in 
models of dynamical symmetry breaking [58], and in little Higgs models [59]. In 
supersymmetric models, an extra U(l)' can provide an elegant solution to the [i 
problem [60j[61], with an effective \x parameter generated by the VEV of the Standard 
Model singlet field S which breaks the U(l)' symmetry. This is somewhat similar to 
the effective \x parameter in the NMSSM [42]. However, with a U(l)' there are no 
extra discrete symmetries or domain wall problems. A closely related feature is that 
the MSSM upper bound of M z on the tree-level mass of the lightest Higgs scalar is 
relaxed, both in models with a U (1)' and in the NMSSM, because of the hSH u Hd term 
in the superpotential and the U(l)' D-term [63]. Light Higgs masses lower than those 
allowed in the MSSM are also allowed experimentally because of possible mixings 
among the Higgs doublets and singlets [611 [65]. F° r specific U(l)' charge assignments 
for the ordinary and exotic fields in some theories one can simultaneously ensure 
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the absence of anomalies; that all fields of the TeV-scale effective theory are chiral, 
avoiding a generalized \i problem; consistency with gauge coupling unification; and 
the absence of dimension-4 proton decay operators [66] . 

There are stringent limits on an extra Z' from direct searches from the Teva- 
tron [H7J and from indirect precision tests at the Z-pole, at LEP 2, and from weak 
neutral current experiments [68]. The constraints depend on the particular Z' cou- 
plings, but in typical models one requires Mz> > (800 — 1000) GeV and the Z — Z' 
mixing angle ctz-z 1 to be smaller than a few xl0~ 3 . (The specific parameters con- 
sidered here yield a Z' mass of around 1 TeV.) Thus, explaining the Z — Z' mass 
hierarchy without fine-tuning is non-trivial. Two of us with J. Erler proposed a su- 
persymmetric U(l)' model with a secluded [/(l)'-breaking sector (sMSSM) in which 
the squark and slepton spectra can mimic those of the MSSM, the EW symmetry 
breaking is driven by relatively large A terms, and a large Z' mass can be generated 
by the VEVs of additional SM singlet fields that are charged under the U(l)' but do 
not directly contribute to the effective \x parameter [OH] . If these fields are only weakly 
coupled to the SM fields, i.e., by U(l)' interactions and possibly soft SUSY-breaking 
terms, then the scale of VEVs in this sector is only weakly linked to the EW scale. 
In particular, we consider the situation in which there is an almost F and D flat 
direction involving these secluded fields, with the flatness lifted by a small Yukawa 
coupling A. For a sufficiently small value for A, the Z 1 mass can be arbitrarily large. 

In this paper, we systematically study EWBG and its correlation with neu- 
tralino CDM in the sMSSM (for a discussion on EWBG in the thin wall regime in the 
sMSSM, also see [TO]). We construct two explicit anomaly-free models. In Model I, 
we embed the SU(3)c x SU(2) L x U{l)y x U(l)' gauge symmetry into a subgroup of 
E 6 . We consider three families of the SM fermions as arising from three fundamental 
27 representations of Eq. The 27's also contain candidates for the Higgs doublets, 
SM singlets (needed to break the U(l)'), and additional exotic fields. Additional 
vector-like pairs from 27 + 27* representations can also exist at low energy without 
introducing anomalies. Consistency with gauge coupling unification is achieved by 
introducing one pair of vector-like Higgs-type doublets from 27 + 27* |7lffi while 
the necessary U(l)' breaking and the generation of masses for the exotic fermions 
is achieved by introducing three pairs of vector-like SM singlets from two pairs of 
27 + 27*. The spectrum of exotics is rather complicated in Model I. In contrast, in 
Model II we add the minimal number of exotic particles with rational U(l)' charges 
to cancel the U(l)' anomalies. This model has the minimal exotic particle content 
but does not respect the simplest form of gauge coupling unification. We present the 
general superpotential and the supersymmetry breaking soft terms for both models. 

For concreteness, we will work in Model I (the Eq embedding). First, we discuss 
the one-loop effective potential at finite temperature in the 't Hooft-Landau gauge in 
the MS'-scheme. We calculate the gauge boson, Higgs, neutralino, chargino, squark, 
and slepton mass matrices, and those for the scalar components of the exotic chiral 

2 Thcse are in addition to those in the three 27's. It is not important for our purposes which ones 
correspond to the two MSSM Higgs doublets. 
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superfields. We also calculate the temperature dependent masses for the Higgs fields, 
the longitudinal components of the gauge bosons, squarks, sleptons, and the scalar 
components of the exotic chiral superfields. There are two relevant phase transitions 
as the temperature decreases: the U(l)' symmetry is broken at about 1 TeV, and 
the EW symmetry at the weak scale. We show that there exists strong enough first 
order EWPT because of the large trilinear term AhhSH u Hd in the tree-level Higgs 
potential. Thus, unlike the MSSM, the stop masses can be very heavy compared to 
the top quark mass. The EWPT features in model II or the other possible anomaly- 
free embeddings are similar, because the exotic particles' contributions to the one-loop 
effective potential at finite temperature are suppressed by their heavy masses. 

Subsequently, we discuss possible CP violation introduced by the extended 
Higgs sector of the sMSSM. As first pointed out in [70], unlike the MSSM in which 
there is no CP- violation in the Higgs sector at tree level, both explicit and spon- 
taneous violation (denote them as ECPV and SCPV) can occur. There are five 
complex parameters in the supersymmetry breaking soft terms and only four gauge- 
independent Higgs phase degrees of freedom, which implies that one of these complex 
phases cannot be removed by field redefinition. The CP-symmetry of this tree- level 
potential can therefore be explicitly broken. In addition, as a result of balancing dif- 
ferent terms in the minimization of the tree-level neutral Higgs potential, the phases 
of the Higgs fields can obtain non-trivial VEVs, which implies that CP symmetry 
can also be spontaneously broken. With loop corrections included, there may coexist 
vacua at finite temperature with broken and unbroken EW symmetry. The values of 
the spontaneous CP phases are usually different in these vacua. 

In this work SCPV provides a direct source for baryogenesis, while its mag- 
nitude is mediated by an explicit phase from the secluded sectoiH. These new CP 
sources do not introduce significant new contributions to the electric dipole moments 
(EDMs) of the electron and neutron. After proper field redefinitions, the CP- violation 
phases will only appear in the Higgs mass matrix through soft masses associated with 
singlet components. Numerical estimates show that, for typical parameter values, 
their contributions to EDMs will be about six or seven orders of magnitude smaller 
than the experimental upper limits. These contributions disappear completely in the 
limit with a trivial explicit CP phase, where the spontaneous CP phases are absent 
in the true vacuum but not where EWBG occurs. 

In the early Universe the first order EWPT is realized by nucleating bubbles 
of the broken phase. The dynamical properties of these bubbles, such as the wall 
profile and expansion velocity, can have important influences on the production of 
the baryon asymmetry. We also study their physics in detail. The VEVs of the 
Higgs fields (including both their magnitudes and phases) are space-dependent, as 
one crosses the bubble wall. The bubble wall thickness is estimated by minimizing 

3 We turn off all CP- violating sources beyond the Higgs sector, e.g., the ones from soft gaugino 
and sfermion masses, which are generally used for EWBG in supcrsymmetric models. Though they 
may provide non-trivial contributions to EWBG, these CP phases usually suffer strong constraints 
from electric dipole moments. 
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the action under the kink ansatz (721 [731 EH- Numerical results show that it varies 
from 3 T^ 1 to 30 T~ l as an approximately monotonically increasing function of the 
phase changes of the Higgs fields. In addition, we argue that the wall velocity in the 
sMSSM cannot be larger than that in the MSSM under the same phase transition 
condition and thus should be non-relativistic. This fact implies that EW sphaleron 
processes have more time to occur and hence will enhance the final baryon asymmetry. 

We then systematically study non-local EWBG in both the thin wall and thick 
wall regimes. Non-local baryogenesis has a great advantage over the local one |75j : 
the baryon number violation processes take place in an effective volume extending 
from the wall surface to a region in the symmetric phase, and thus we may expect an 
enhancement in the generation of the baryon asymmetry. According to the bubble 
wall profile, we calculate the contribution from r leptons in the thin wall regime, and 
the ones from stop, charginos and top quarks in the thick wall regime. We find that 
the CP- violating currents induced by stops and charginos are very different from those 
obtained in the MSSM: (1) the current of the Vidv 2 —V2dvi type (in our formalism, the 
concrete expression is a little different) never requires d tan j3 ^ in the bubble wall 
due to the variance of the relevant CP phases crossing the walfl; (2) there is a new 
type of CP- violating current at the leading order which is proportional to v(Oi with 
i — 1,2 and is absent in the MSSM. The CP-violating current of the new type has 
important influence on the EWBG. First, the stop contribution can be quadratically 
enhanced by a large soft Ah t parameter. Second, in addition to r leptons, top squarks 
and charginos, top quarks can also play a significant role in the EWBG. All of these 
features are results of the SCPV-driven EWBG, so they are not sensitive to the 
concrete embeddings of the sMSSM. Numerical calculations show that the produced 
baryon asymmetry is large enough to explain the cosmological observation today. 
Though the EWBG is directly driven by SCPV, there is no dilution problem between 
matter and anti-matter bubbles. Such a dilution is caused by a Z 2 discrete symmetry 
which usually exists in the SCPV scenarios. In our work, however, this symmetry has 
been explicitly broken at tree level by the explicit CP- violating phase. 

We also study the correlation between EWBG and cold dark matter (CDM) 
in the sMSSM. Large trilinear soft parameters (used for strong enough first order 
EWPT) more often arise from gravity-mediated SUSY breaking, where the lightest 
super symmetric particle (LSP) is usually the lightest neutralino. We study the possi- 
bility of describing the EWBG and neutralino CDM in the same framework. Though 
we do not scan the whole parameter space, we find that there indeed exist some re- 
gions where strong enough first order EWPT, large CP phase variations across the 
bubble wall, reasonable baryon asymmetry, as well as acceptable neutralino LSP relic 
density can be achieved simultaneously. The relevant neutralino mass m x o is close to 

4 In the MSSM, there also exists a CP-violating current of the V\dv2 + V2dv\ type for charginos 
which arises from a resummation of corrections associated with higher order Higgs insertions or 
multiple scattering effects in the bubble wall. Though not suppressed by 9tan/3, the current of this 
type is usually subdominant unless d tan/? is suppressed (in our model it is always subdominant). 
For relevant or more general discussions, see [751 IZ3 HS1 HZ1 [Ml 1301 [2D 132] ■ We will therefore not 
discuss the currents of this type. 
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Mz/2. Finally, we comment briefly on possible cosmological signals: superconduct- 
ing cosmic strings and gravitational waves (GWs). We expect that particle emission 
from the decays of cosmic strings and GWs from EWPT may be observed within the 
foreseeable future. 

This paper is organized as follows. We briefly review the tree-level Higgs 
potential in the sMSSM and construct two anomaly-free embeddings in Section 2. 
The one-loop effective potential at finite temperature is considered in Section 3. In 
Section 4, we show that the EWPT is strongly first order. We study CP violation in 
Section 5 and the bubble wall physics in Section 6. In Section 7, we systematically 
study non-local EWBG in both the thin wall and thick wall regimes, and discuss the 
contributions from leptons, squarks, charginos and quarks in detail. In Section 8, we 
illustrate that there exists common parameter space where the baryon asymmetry 
and the CDM can both be explained by EWBG and neutralino LSP, respectively. 
Some simple comments on the cosmological signals of cosmic string decays and GWs 
are given in Section 9. The last section is our discussion and conclusions. 



2 Two Anomaly- free U(l) f Models 

Let us briefly review the sMSSM [69]. There are one pair of Higgs doublets, H u and 
Hd, and four SM singlets, S, Si, S2, and S3. The U(l)' charges for the Higgs fields 
satisfy 

Q' Hd + Q' Hu + Q' s = 0, Q' s = -Q' Si = -Q' S2 = (6) 

The superpotential for the Higgs is 

W H = hSH u H d + XSiS 2 S 3 , (7) 

where the Yukawa couplings h and A are respectively associated with the effective fi 
term and with an (almost) F and D-flat direction. For simplicity, we assume that 
terms such as SfS^, Sf 1S3 and their associated soft terms are absent, and that there 
are no bilinears SSi )2 in Wh- The existence of a number of SM singlets and the 
non-diagonal nature of the superpotential is in part motivated by explicit superstring 
constructions [78]. The corresponding F-term scalar potential is 

v F = /. 2 (i^ o i 2 ii/ n o r+i^i 2 i^T+i^i 2 i^ o i 2 ) 

+\\\Si\ 2 \S 2 \ 2 + \S 2 \ 2 \S,\ 2 + \S,\ 2 \Si\ 2 ) . (8) 
The D-term scalar potential for the neutral fields is 

Vd = ^(I^°| 2 -I^°| 2 ) 2 

+\g% (^l^l 2 + ^J^T + Qkl^T + E^I^I 2 ) > W 



s 



where G 2 = g\ + g\\ gi,g2, and gz> are the coupling constants for U(l), SU(2)l and 
t/(l)'; and Q' is the U{1)' charge of the field <p. 

In addition, one introduces the supersymmetry breaking soft terms 

3 

mH 2 I tjO 1 2 i 2 I ttO 1 2 , 2 I cl2 , 2 | c 1 2 

Ko/t = m H d \H d \ + m H JH u \ +m s \S\ + }^ m spi\ 

i=i 

- (A h hSH° d H° u + A X XS 1 3 2 S 3 + m 2 SSi SS 1 

+m 2 SS2 SS 2 + m 2 SiS2 S\S 2 + H.C.) . (10) 

rrigg., i — 1, 2 are needed to break two unwanted global t/(l) symmetries (for a recent 
discussion of such global symmetries, see [79]). There is an almost F and D flat 
direction involving Si, with the flatness lifted by a small Yukawa coupling A. For 
a sufficiently small value of A, the Z' mass can be arbitrarily large. For example, 
if h ~ 10 A, one can generate a, Z — Z' mass hierarchy in which the Z' mass is of 
order 1 TeV. The dimensional parameters in V s ^ t are specified in arbitrary units, 
and then rescaled after the (effective) potential is minimized so that the EW scale is 
v = y/vl + v 2 ~ 174 GeV, with v x = (H§ and v 2 = (H$). 

In (69], m 2 s s was set to for simplicity. In that case, the fields can be defined 
so that Ah, A\, rn 2 SSv and m 2 SS2 are all real and positive. The minimum then occurs 
at a point at which all of the fields are real and positive, i.e., there is no CP violation 
in the tree-level Higgs potential. For m| S2 ^ 0, however, there will in general be CP 
violation. 

The U(l)' charges for the fermions and the anomaly cancellations were not 
discussed in Ref. [69], so in the following two subsections we construct two anomaly- 
free models, one of which is consistent with minimal gauge coupling unification. 

2.1 Model I: sMSSM with E 6 Embedding 

U(l)' models necessarily imply that new fermions are needed for anomaly cancella- 
tions [55l [66]. Since all complete representations of E§ are anomaly-free, it is con- 
venient and conventional to simply consider the U(l)' charge assignment and exotic 
particle content of Eq as an example of an anomaly-free construction. We are not 
considering a full E e grand unified theory, because a light Z' would prevent a large 
doublet-triplet splitting, leading to rapid proton decay if the E 6 Yukawa relations 
were enforced. (Detailed studies of E 6 theories with broken Yukawa relations may be 
found in jgQ|.) 

The Eq gauge group can be broken as [HH E2] 

E 6 -> SO(10)x [/(l)^ -> SU(5)x U(l) x x U(l)^ . (11) 

The U(l)ip and U(l) x charges for the Eq fundamental 27 representation are given in 
Table [TJ The U(l)' in Model I is one linear combination of the U(l) x and U(l)^ 

Q' = cos9 Q x + sm9 , (12) 
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Table 1: Decomposition of the Eq fundamental 27 representation under SO(10), 
SU(5), and the U(l) x , U{1)^ and U(l)' charges. 



50(10) 


SU{5) 


2VWQ X 


2^^ 


2^0, 


16 


10 {u, d, u, e) 


-1 


1 


-1/2 




5 (d, v, e) 


3 


1 


4 




IN 


-5 


1 


-5 


10 


5 (D, H' u ) 


2 


-2 


1 




5(D,H' d ) 


-2 


-2 


-7/2 


1 


1 5 L 





4 


5/2 



where the angle # will be chosen to ensure the needed U(l)' charges for the SM 
singlets. For simplicity, we assume that the other U(l) gauge symmetry from the 
orthogonal linear combination of the U(l) x and U{1)^ is absent or broken at a high 
scale. 

We assume three 27s, which include three families of the SM fermions, one 
pair of Higgs doublets (H u and H d ), and a number of SM singlets, extra Higgs-like 
doublets, and other exotics. The embeddings of the SM fermions are obvious. For 
defmiteness, we assume that the Higgs doublets (H u and Hd) are the doublets in 10 
(or 5 and 5) in the third 27. (Of course, the exact identification of the H u , Hd and S 
with specific 27s is irrelevant.) In addition, we assume that the four SM singlets S, 
Si, S2, S 3 are the Sl, SI, SI and N*, respectively, in two (partial) pairs of 27 and 
27*. To cancel the U(l)' anomalies, we introduce X and X 3 , which are Sl and N, 
respectively, in the two pairs. Thus, (S, Si), (X, S 2 ) and (X 3 , S 3 ) are three pairs of 
vector-like SM singlets and will not introduce any anomalies. To be consistent with 
minimal gauge coupling unification, we introduce one pair of vector-like doublets H' u 
and H' u from a pair of 27 + 27* I?!]. For simplicity, we assume that the other particles 
in the 27 + 27* pairs are absent or are very heavy and decouple at low energy. From 
the requirement Qs = \Qs^ Qs L = \Qn*i we obtain 



tan # = 1^ . (13) 
9 V ' 

The exotic particles in Model I are Di, Di, S L , Ni, where i = 1,2,3; H£ and H' d k , 
where k = 1,2; X, X 3 ; and H' u and H' u . The U(l)' charges for the Standard Model 
fermions and exotic particles are also given in Table [TJ 
The general superpotential is Wh + Wy, where 

Wy = aSXX 3 + agSSiNj + rfKK + ^.O.H.a;, + h%H d Qid 3 + h%H d U^ 
+h^H u L l N ] + apD t D 3 + a%SH' u k H' d l + X^Q.Ku, 
+\? j H' u L l N j + \) jli: I),a j< ,„ + % jm DiQjL m + \% m Dd 3 N m 
+X A ijm D l Q J Q m + X\ jm DiUjd m + X e ki 4>S k S l L + M^ 2 , (14) 



10 



with Qi = (v,i,di), Li = (b>i,ei), and i,j,m = 1,2,3 while k, I = 1,2. For simplicity, 
we will assume that the Yukawa couplings h%, hfj, h\- and hf- and those associated 
with the exotics are diagonal in the numerical calculations. In addition, to explain 
the neutrino masses and mixings, we introduce a SM singlet which is neutral under 
U(l)'. Integrating out <fi, we obtain the non-renormalizable terms SkSiS^S^/M^, and 
then Majorana mass terms for S^S^ will appear after U(l)' breaking. Thus, the 
light neutrino masses and mixings in Model I can be generated via a double-see-saw 
mechanism naturally [83| [HI] . We see from Eq. Q14p that all of the exotics except 
H' u and H' u can have masses generated by the VEV of S. The latter can be given a 
super symmetric mass \J . This reintroduces a form of the fi problem, though not for 
the Higgs fields associated with EW symmetry breaking, and is clearly a flaw of the 
specific construction. However, this segment of the model is only loosely connected 
with the issue of concern in this paper, i.e., EWBG, so we will tolerate it rather than 
going to a more complicated model. 

The new scalar supersymmetry breaking soft terms (in addition to Eq. (fTUl) ) 

are 

V s A oft = -A^SXXs-A^apSiNj-B^^H^-A^Q^uj 

-A^hfjHjQidj - Ahe.h^HdLiEj - A h Nh^ j H u L l N j - A^apD^ 

- A ^ m ^jm D iQjQm - A| m A%4 - A^Xt^SkSl 

-B^<f) 2 + H.C., (15) 

3 
i=l 

+m%]R i \ 2 + m\ \S^ + m% | A| 2 + ™|JA| 2 ) 

k=l 

+m||X| 2 + m 2 ?3 |X 3 | 2 + m^|<| 2 + m|,|^| 2 , (16) 

where we have assumed that the mass-squared terms are diagonal. We assume that 
only S, Si, S2, S3, H®, and H% acquire VEVs. 

The vector-like particles H' u and H' u can decay via the Yukawa Xf, and \fj 
terms in Eq. (j!4p . In addition, the charged exotic particles Di and Dj can decay 
through the Yukawa A" m terms with n = 1, 2, 5 in Eq. (fT4"j) [521 [55] . However, 
to avoid the proton decays via dimension-5 operators, we must forbid some of the 
Yukawa A™ m terms. For example, we can choose A™ m = either for n — 1, 2, 3 or 
for n — 4, 5. 
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2.2 Model II: Minimal Anomaly-free Model 

In this subsection, we consider the minimal sMSSM in which all the particles have 
rational U(l)' charges. In order to cancel the gauge and gauge-gravity mixed anomaly 
due to the U(l)' gauge symmetry, we introduce the exotic particles D\, D\, D 2 , 
D 2 , V, V, Xi, Xi, X 2 , X 2 , X 3 , X 3 , Xi, and X 4 . The quantum numbers under 
SU(3) C x SU(2) L x U(1) Y x U(l)' for the SM fermions, Higgs fields and extra exotic 
particles are given in Table [2j If one allowed irrational U(l)' charges, only D\, D 2 , 
D 2 ,V, and V would be needed. 

With the U(l)' charge assignments for the particles in Table [2j the conditions 
for [SU(3)c] 2 U(l)' and [Gravity] 2 £7(1)' are automatically satisfied, and the condition 
for [SU(2)l] 2 U(1)' is the same as that for [U(1)y} 2 U(1)' . The anomaly-free conditions 
for [SU(2) L } 2 U(1)', U(l) Y [U(l)'} 2 and [U(l)f, respectively, are 

Cl + c 2 + 3(a + 36) = , (17) 

-3a 2 + 36 2 + 3(a + ci) 2 - 6(6 + c 2 ) 2 + 3(b + ci) 2 -c\ + c\ (18) 
-d\ + (s + d x f -d 2 2 + (2s + d 2 f -dj + (2s + d 3 ) 2 = , 

6a 3 - 3(a + c 2 f - 3(a + ci) 3 + 186 3 - 9(6 + c 2 ) 3 - 9(6 + ci) 3 (19) 
+3df - 3(s + di) 3 + 3d 3 , - 3(2s + d 2 ) 3 + d 3 - (2s + d 3 ) 3 
+2d 3 - (d 4 - s) 3 - (d 4 + s) 3 + 7s 3 + 2c 3 + 2c 3 = . 

The simplest solution is c\ = 9/10, c 2 = 9/10, s = —9/5, a = —18/5, b — 1, d\ — —2, 
d 2 = -1/5, d 3 = 9/5, d 4 = 27/10, d 5 = 9/10, and d 6 = -3. 
The general superpotential in Model II is 

W = a 1 SD 1 D 1 + a 2 S 3 D 2 D 2 + a 3 S 3 VV + a 4 5 1 X 1 X 1 + a 5 S 2 X l X l 
+a 6 SX 2 X 2 + a 7 S , iX 3 X3 + a$S 2 X 3 X 3 + (X9SX4X4 
+a[SiX±Xi + a 2 5'iXiXi + ^^XiXx + ck^^-^i-^i 
+h u tJ Q i H u u J + hfjHaQidj + h^HaLiEj + hf j H u L i N j 
+\} j D 1 Q i Q J + Xijhuid, + A^Xl + A^e.Xi . (20) 

All of the exotics can be made massive by U(l)' breaking. 

The charged exotic particles (D iy Dj), and (V, V) can decay through the 
Yukawa A^-, A 2 ,-, A 3 and A^ terms in Eq. (EDj) [82l [85] . where we have assumed that 
my > mx,. After the U(l)' breaking, the vector-like particles (Xi, Xi) are neutral. 
There is no discrete symmetry after U(l)' breaking, so the vector-like particles (Xi, 
Xi) can decay via higher dimensional operators. 
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Table 2: Minimal Anomaly-free Model. Quantum numbers under SU(3)c x 
SU(2) L x U(1) Y x U(l)' for the left-chiral Standard Model fermions 
(Qi, Ui, di, Li, Ni, e;),_Higgs fields ( H u , H d , S, Si, S 2 , S3), and extra 
exotic particles ( D\, D\, D2, D 2 , V, V, X\, X\, X2 and X2). 



Particles 


Quantum Numbers 


Particles 


Quantum Numbers 


U 


(1; 2; -1/2; a) 




(3; 2; 1/6; b) 


N t 


(1; l;0;-(a + c 2 )) 


Mi 


(3; 1; -2/3; -(b + c 2 )) 


ei 


(1; 1; 1; -(a + ci)) 




(3; 1; 1/3; -(6 + Cl )) 


H d 


(1; 2; -1/2; c x ) 


H u 


(1; 2; 1/2; c 2 ) 


S 


(1; 1; 0; s) 


S3 


(1; 1; 0; 2s) 


Si 


(1; 1; 0; -s) 


S 2 


(1; 1; 0; -s) 


Di 


(3; 1; -1/3; d l ) 


Di 


(3; 1; 1/3; -(s + di)) 


D 2 


(3; 1; -1/3; d 2 ) 


D 2 


(3; 1; 1/3; -(2s + d 2 )) 


V 


(i; i; -i; d 3 ) 


V 


(1; 1; 1; -(d 3 + 2s)) 


Xi 


(l; l; 0; d 4 ) 


Xi 


(1; 1; 0; -(d 4 - s)) 


x 2 


(1; 1; 0; d 5 ) 


x 2 


(1; 1; 0; -(4 + a)) 


x 3 


(1; 1; 0; d 6 ) 


X3 


(1; 1;0; -(d 6 - s)) 


x. 


(1; 1; 0; d 6 ) 


x. 


(1; 1; 0; -(d 6 + s)) 



3 One-Loop Effective Potential at Finite Temper- 
ature 



In the 't Hooft-Landau gauge and in the MS'-scheme, the one-loop effective potential 
at finite temperature is [86J 

V cS (<p,T) = V Q (<f>) + V 1 (<j>,0) + AV 1 ( ( f),T) + AV daisy ( ( f),T) , (21) 

where Vo(4>) is the tree-level potential, Vi(4>, 0) is the one-loop zero temperature cor- 
rection, AVi(0, T) is the temperature dependent one- loop correction, and AVd a is y (</>, T) 
is the multi-loop daisy correction. The explicit expressions for Vo(0), Vi{4>, 0), AVi{4>, T) 
and AV d aisy(0, T) are 

W) = V F + V D + V s H oft , (22) 



(23) 



AT4(0,T) 



rpl 

2^2 



^ 2 (0) 

J^2 



(24) 
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daisy I 



T) 



T 



12tt 



rii [M^ ((f), T) 



(25) 



i=bosons 

where mj(0) are field-dependent masses, rij are the number of degrees of freedom, and 
Ci are constants dependent on the regularization scheme: in the MS-scheme that we 
are assuming, Ci = 5/6 for the gauge bosons and 3/2 for scalars and fermiona^]. The 
function Jj comes from the one-loop corrections to the effective potential at finite 
temperature [86J. For bosons 



J B (m 2 (</>)/T 2 



dx x 2 log 



. L.2 i rn 2 (4>) 



(26) 



At relatively high temperature, < 2.2, one can expand 



J B (m 2 (0)/T 



2\HT 



_7T* 7T 2 

45 + 12 T 2 



m 2 ((j)) 7T ( m 2 (0)\ 



6 



3/2 



T 2 I 



i "A0) log ^_ 27r 7/ 2 g ( 



32 T 4 



a fe T 2 



C(2£+l) ( 1\ f m 2 (0) \ 

(£ + i)! v 2/ v 471-2712 / 



^+2 



(27) 



where a& = 167r 2 exp(3/2 — 27^) (loga^ = 5.4076) and ( is the Riemann (^-function. 
At relatively low temperature, > 2.2, we can expand [8 



J s (m 2 (0)/T 



2\LT 



45 + (2 
15 T 



m(0)\ 



3/2 



-m(0)/T 



1 + 



m((f)) 



The function for fermion ^ with mass m is 



J F (m 2 (V0/T 2 



x 2 log 



x 2 + = iW 



(2? 



(29) 



For ^ < 1.8, one finds [86] 



J F (m 2 (^)/T 



2\.ffT 



360 

7T 7 



1 m^iip) m 2 (ifj) 
log- 



7r 2 m 2 (^) 

24 T 2 32 T 4 a/T 2 



e=i 
1 
2 



m 2 (ip) 

7T 2 T 2 



£+2 



(30) 



5 For the £>i?-scheme, Ci in Eq. (|23|) is 3/2 for gauge bosons, scalars, and fermions. Our discussions 
and conclusions do not depend on the scheme. 
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where a/ = it 2 exp(3/2 — 2^e) (logo/ = 2.6351). At low temperature, > 1.8 [87] . 



M^m/rr = ^ + (i) 1/2 (^) 3/2e """* )/T 

15 T 
1 H r~- + ... 

8 m(^) 



(31) 



The m(</>) and ra(^>) are the Higgs field dependent masses from the eigenvalues 
of the mass matrices. Thus, we need the mass matrices for all the particles in the 
models. The tree-level mass matrices^] for the particles in Model I are presented in 
Appendix A. (Those for the ordinary and exotic fermions are trivial.) 

There are small discontinuities between the Js(m 2 (0)/T 2 ) ffT and Js(m 2 (0)/T 2 ) LT 
at ~ 2.2, and between the J F (m 2 (0) /T 2 ) HT and J F (m 2 (0)/T 2 ) LT at ^ ~ 1.8. 
When we numerically calculate the minimum of the potential, we need smooth inter- 
polations for JB(m 2 ((j))/T 2 ) and Jj?(m 2 (0)/T 2 ). Thus, in our numerical calculation 
we use the following approximation 

t ^r-tf)ti x MmH4>)/T > r (32) 



Mm\m 2 ) = taDh (- m ?w/ r2 + 1.8') + 1 x Jr(m 2 W/ry T 

+ ^r-i8i + i x jF(m , W)/rV T (33) 

The temperature dependent scalar mass-squared A4 2 (4>,T) is obtained from 
the m 2 (<j)) by adding the leading temperature dependent self-energy contributions 

n*(r), 

X 2 (0,T) = m 2 (0) + n <I> (T) , (34) 

where 11$ (T) is proportional to T 2 [SS]. For the gauge bosons, only the longitudinal 
components receive such contributions. The n$(T)'s for the particles in Model I are 
presented in Appendix B. The daisy corrections are of higher order, but important 
at high temperature. 



6 The use of tree-level mass matrices leads to a slight inconsistency when they are evaluated at the 
minimum of the full one-loop effective potential. In particular, the tree-level approximations to the 
Goldstone Boson masses are sometimes shifted slightly from zero. We avoid singularities by setting 
slightly negative mass-squares to zero. We have checked that these small shifts have negligible effects 
on the phase structure, as should be apparent from the formula (|27|) relevant for small m lf . 
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Table 3: Three sets of typical parameter values. The energy units are 69 GeV, 92 
GeV and 120 GeV, corresponding to m? Q = 5.0, 9.0 and 13.0, respec- 
tively. In this paper, we will denote them b and c. Note that 

the soft mass squares in the superpotential are relatively small, i.e., the 
symmetry breaking is ^4-term dominant. M[, Mi and M 2 denote the 
three neutral gaugino soft masses. 



h 


A h 


A 


A x 


m S Sl 


m ss 2 






0.8 


4.2 


0.06 


3.3 


0.02 


0.1 


4.8 xl0~ 4 


5.0 / 9.0 / 13.0 


m 2 H o 


m s 


m k 


m s 2 


m| 3 


M( 


Mi 


M 2 


-0.1 


0.5 


0.03 


0.03 


0.03 


1.5 


1.5 


3.0 


m ' 2 n 
Q3 


2 

m- 3 


m '% 
4 


Ql,2 


% 


d R 






8 


8 


25 


25 


25 


25 


25 


25 



4 The Electroweak Phase Transition 

Generally, the EWPT in the sMSSM is more complicated than those in the SM, 
MSSM, and NMSSM due to the secluded L r (l)'-breaking sector. The temperature de- 
pendent corrections to the effective potential are a function of the masses Aii((f>,T). 
At low temperature, the temperature dependent corrections are negligible, while at 
very high temperature they are dominant. At the critical temperature, the temper- 
ature dependent corrections are comparable to the effective potential at zero tem- 
perature. Our one-loop effective potential is a function of 10 physical variables (the 
complex VEVs of the six neutral Higgs field^ minus two gauge degrees of freedom) . 
There might exist local minima near the critical temperature. In our analysis we 
always choose the global minimum (at one-loop level plus daisy correction). 

Because the Z' mass is assumed to be of order 1 TeV, the VEVs of the SM 
singlets Si, 5*2, and S3 are about an order of magnitude larger than those of the Higgs 
doublets and S [69]. There are therefore two phase transitions: the U(l)' symmetry 
breaking at TeV scale and the EW symmetry breaking at the weak scale. 

In Figure [lj we plot the Higgs VEVs versus temperature for the given param- 
eters values in case b (see Table EJ). From the figure, one can see that the U(l)' phase 
transition occurs just below 1200 GeV, and the EWPT at about 120 GeV. We are 
interested in whether the EWPT is strongly first order. From the figure, we also see 
that the VEVs of H%, H® and S all experience a first order transition due to their 
discontinuity at the weak scale. The magnitudes of the VEVs of the Higgs fields are 
plotted at the global minimum of the potential for the given temperature. There is an 
abrupt change in the Si VEVs at temperatures around 400-450 GeV because another 

7 We ignore the possibility of charge or color breaking minima. 
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Figure 1: The Higgs VEVs vs. temperature T. The VEVs of H^ d and S are 
multiplied by 10. The VEVs of Si, S2 and S3 are indistinguishable. 

local minimum becomes the global oneEl 

After the U(l)' symmetry is broken, Si, S2, and S3 have non-zero VEVs of 
order 2 TeV, while the VEVs for H% and H® are still zero. The supersymmetry 
breaking soft terms m 2 SSi SSi and m 2 SS2 SS2 generate a linear term for S in the po- 
tential after the U(l)' symmetry breaking, which induces a non-zero VEV for S. At 
high temperature, (S) is tiny because m 2 SSi and m 2 SS2 are small and the temperature 
dependent corrections to the S mass are large. At lower temperature, (S) becomes 
larger, as can be seen from Figure [D At the EWPT the term AhhSH u Hd induces a 
first order phase transition in the VEV of S as well. 

The first order EWPT can be strong because of the trilinear term AhhSH u Hd 
in the tree-level potential. Unlike in the MSSM [301 ED E2], in which the phase 
transition is associated with loop effects, there is no upper bound on the lightest stop 
mass. This feature is independent of the concrete embeddings of the sMSSM, though 
we are using an effective potential of finite temperature in Model I as an illustration. 
This feature is also shared by the NMSSM H HQ M> HE] , the nMSSM [EH [52], the 
U(l)' models with no secluded sector [89], and the singlet extensions of the SM [2~T] . 

In Table IU we calculate v(T c )/T c versus the stop soft mass square m ~3 
and Ah while the other parameters are fixed to the values given by Table |3j with the 

8 In a real transition the fields would most likely remain in the same minimum, but that would 
have little effect on our conclusions. Detailed investigations of issues such as these are very model 
dependent and beyond the scope of this investigation. 
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Table 4: v(T c )/T c versus the stop soft mass squares mj~ , m~ 3 and Ah, where the 

other parameters are fixed to the values in case b (see Table [3]) . For 
simplicity, we assume m% = to? 3 . The scales are different for different 

Q'A U R 

points, but typically they are about 100 GeV. The slashes in the table 
mean that there is no first order EWPT. 





2.0 


3.0 


4.0 


5.0 


4.0 


/ 


1.29 


1.59 


2.06 


8.0 


0.91 


1.02 


1.21 


1.50 


12.0 


/ 


0.56 


1.03 


1.22 


16.0 


/ 


/ 


/ 


1.37 



choice m 2 H0 = 9.0. The scales are different for different points, but typically they are 

d 

about 100 GeV. As expected, when Ah decreases (with fixed mj~ 3 and m? 3 ), v(T c )/T c 
decreases. The first order EWPT is not strong enough when Ah < 2 — 3. This table 
also displays the dependence of v(T c ) /T c on mj~ g and m? 3 . For fixed A h , the strength 

of the first order EWPT usually becomes weaker as mj~ increases. Unlike the MSSM, 

even for heavy stop soft masses (or physical masses) one can still have a strong enough 
first order EWPT. 

For simplicity, we did not turn on any CP phases for the discussion in this 
section. We will see later that the results obtained in this section are insensitive to 
these phases. 



5 Generation of CP Violation 



5.1 CP Violation at Tree Level 

One nice feature of the sMSSM is that both ECPV and SCPV can happen in the 
Higgs sector at tree level [70]. This model involves six complex neutral fields. To 
study the CP- violation, we define 

Hi = \H° d \e^ , H° u ee \H° u \e^ , S = \S\e* a , 

S x = , S 2 = \S 2 \e ia2 , ^ = l^le^ 3 . (35) 

Among these six phase degrees of freedom, only four of them, which are given by 

(3\ — a + ai , p 2 = « + «2, 
f3 3 = a + 9i + 9 2 , At = «i + a 2 + a 3 . (36) 
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are invariant under the U(l)y and U(l)' gauge transformations. The two remaining 
degrees of freedom 



Z'iYfcOfc , HiQfcdfa , (37) 

correspond to Goldstone bosons associated with the two spontaneously broken gauge 
symmetries, and decouple from the Higgs sector. Here, and Q'± are the U(l)y and 
U(l)' charges for fa, 9^ is the phase of fa, and the index 11 fa" runs over all six complex 
Higgs fields. The complex phases of the Higgs fields 9^ are linear combinations of the 
gauge-invariant and gauge- dependent components. Of course, physical consequences 
only depend on the gauge-invariant components. 

There are five supersymmetry breaking complex parameters in the tree-level 
Higgs potential Eq (IBl TTOl) . Ayjn, A\\, m 2 SSi , m 2 SS2 and m| s . Four of the five complex 
phases can be resolved by 2,3,4- Without loss of generality, we assume that 

m\ S2 = KsJe^. (38) 

is complex and the other parameters are real and positive. If 7 is not equal to zero 
or 7r, there will be ECPV at the tree level. Unlike the MSSM, SCPV can also occur 
at tree level in this model. Let us rewrite the supersymmetry breaking soft terms for 
the neutral Higgs bosons as 

K H oft = ^o|^| 2 + m^o|^| 2 + m||5| 2 + X:m|j^| 2 

i=l 

-2A h h\S\\H° d \\H°\ cos/3 3 - 2A X X\S 1 \\S 2 \\S 3 \ cos/5 4 

-2m| Sl |5'||5'i| cos/5i - 2m 2 SS2 \S\\S 2 \ cos/3 2 

-2 \m 2 SlS2 1 1 S ± I \S 2 1 cos(- + fa + 7) ■ (39) 



It is easy to see that a mildly dominant Tn 2 SlS soft term will make (3\ or (5 2 7^ 0, 7r for 
7 7^ 0, 7r and hence lead to SCPV. Meanwhile, the strength of SCPV can be controlled 
by the explicit CP phase 7. 



5.2 CP Violation at Finite Temperature 

Explicit CP phases from soft parameters beyond the Higgs sector could also con- 
tribute to EWBG. However, it is advantageous to turn off these phases and use the 
spontaneous CP phases as the source for baryogenesis. To implement the EWPT, we 
need to include radiative corrections to the Higgs potential, as discussed in Section |3j 
At finite temperature there may coexist multiple vacua, in which spontaneous CP 
phases usually have different values. If the phase transition between two vacua is 
realized through bubble nucleation, as in a first order EWPT, these CP phases are 
space-dependent as one crosses the bubble wall. In the next sections we will see that 
the generated baryon asymmetry depends on these variations. 
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To make this more explicit, it is convenient to write the particle interactions 
during EWBG using the individual Higgs phases. We define the two gauge- dependent 
phases as A and B, given by 



Yh°0i + Y H °j2 + Y s a + EiYs.on = A , 
Q'hA + Q'hJi + Q's* + ^Q's^ = B ■ ( 4 °) 
Inverting Eq. ( l36l) and Eq. (140|) . we obtain 

ft = -^i-iA + lft + ^-lix-As, 

117 2 14 4 

2 3 1 4 3 8 

„ 2 = __ /n + _ A __ A + _ A + _ i4 __ B , 

°» = >4 ft+ ^ + ^-^ + I B - (41) 

where we have used the rescaled £7(1) y and C/(l)' charges for the Higgs particles in 
Model I 

Y H o = -Y H o = -l,Y s = Y Sl = Y S2 = Y S3 = 0; (42) 

Q'u d = ^ , Qk = 75 > Q's = -Qk = -Qk = \q's 3 = \ • (43) 

Because A and B correspond to the Goldstone bosons of the U(1)y and U(l)' break- 
ing, they can be rotated away by gauge transformations. 

The variations of the Higgs phases crossing the bubble wall are given by 

A<**> = - (0*),, (44) 

where are the values of the Higgs phases in the false and true vacua, respec- 

tively. More explicitly, 

A0! = A6 2 = - l -Af3 l - l -Af3 2 + ?-Af3z + ^-A{3 i 
5 5 15 15 

Aa = ?A/?1 + \a{S 2 + - -^A/3 4 , 

5 5 15 15 

Ac*i = ?A/?1 - ^A(3 2 - 1a/3 3 + ^A/? 4 , 
5 5 15 15 

Aa 2 = —API + ^A/5 2 - iAft + ^-A(3 A , 
5 5 15 15 

A« 3 = -^A/51-iA/5 2 + ^A/53 + ^A/? 4 . (45) 
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These quantities are not uniquely defined, due to the uncertainties in the winding 
numbers of the (Pi) ft- If there is a large difference between the winding numbers of 
(Pi) f and then AO^ can be large, and hence the variation of 6^ in the bubble 

wall is also large. However, the nucleation rate for a phase transition between two 
such vacua is generically small [70], so we will require \Pi\ < n. 

Let us consider how AO^ is affected by the behavior of the related field vari- 
ables in the wall. For a qualitative discussion, we neglect all loop as well as finite 
temperature corrections, even though they are important for the complete physical 
picture. Due to h 3> A, the \Si\ are fixed to almost the same large VEVs in both the 
symmetric and broken phases, which leads to (P±) f )t = and A/? 4 = 0. The change in 
P 3 , however, is subtle. For the unbroken phase, V\ = v 2 = and hence P3 is not fixed. 
But for non-local EWBG the relevant scattering processes mainly occur in the bubble 
wall where f 12 do not vanish, so we can understand j3% in the unbroken phase as its 
asymptotic value as v\ and v 2 approach zero. From the tree-level potential Eq. f [39l . 
it is easy to see that the non-vanishing Higgs VEVs lead to /3 3 = inside the bubble 
wall except the boundary region between the wall and the false vacuum. We simply 
assume A/?3 = 0. The phase changes of the Higgs fields therefore are mainly due to 
A/?i and Aj3 2 - To illustrate how to obtain large APi and A/5 2 , let us consider their 
first order differential equations from the tree- level potential Eq. fl39|) . 

(m 2 SSl \S\\Si\ cos/3i + |m| i52 ||S , i||S , 2| cos(-/3i + p 2 + 7)) A/?i 
-I'^liSall'S'ill^l cos(-/?i + p 2 + t)A/3 2 = -m| Sl |Si| sin / 9iA|5'| 

(m 2 SS2 \S\\S 2 \ cos/3 2 + Iml^JlSiHSy cos(-/3i + p 2 + 7)) Ap 2 
-\m 2 SlS2 \\Si\\S 2 \ cos(-/3 1 +/3 2 + 7 )A/? 1 = -m 2 SS2 \S 2 \sinp 2 A\S\. (46) 

Obviously, large AO^ or large A(3\ and Aj3 2 requires large A\S\, and also depends 
significantly on 7. 

Let us consider a parameter subspace: lOOlms^ | 2 > m 2 SSl ~ 0-l m 55 2 > 
1 1 ?Ti. ,s 2 1 2 - Assume that (S) / (Si t2 ^) is small (~ 0.01) in the false vacuum and large 
(~ 0.1) in the true vacuum. In the broken phase, the m 2 SSi and m 2 SS2 terms will 
dominate over the m 2 SiS2 term in the tree-level potential and thus favor small (3\ and 
f3 2 . In the symmetric phase, the dominant terms change to m 2 SS2 and rn 2 SiS2 due to 
the sharp decrease of the VEV of S. For 7 ~ ir, Pi is pushed to a large value in 
this region. A large A/3i is then generated. We will follow this line to generate large 
AOfa in our analyses. In Figure El we illustrate the 7-dependence of Pi and p 2 , in the 
false and true vacua respectively. The relevant 7-dependence of A6^. is illustrated in 
Figure [31 For both of them, the parameter values are from case b (see Table [3]) . 

Besides the baryon asymmetry, it would be desirable to also account for the 
CDM. In the scenario of the lightest neutralino CDM candidate, the CDM 

generation has a strong dependence on the mass (m x o) and composition of Xv Since 
m x o is sensitive to the soft mass terms of H®, H® and S due to their roles in deter- 
mining tan @, we will consider EWBG with respect to 7 in three cases first ( cases a, 
b and c) which are specified by different m 2 H0 values (see Table [3]), and then discuss 
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Figure 4: A^i^, A a vs. 7(vr), in cases a, b and c. 

the continuous dependence of EWBG and the Xi re li c density on m 2 H o for a given 7 
value. In Figure H] we show the 7 dependence of A6x = A6 2 and Aa in cases a, b 
and c (since the baryon asymmetry generation mainly couples to these three phases). 
From this figure we see that for a given 7 value a larger m 2 value gives a smaller 

phase variation crossing the wall. This is because a large m 2 H0 value make \S\ small 

d 

in both vacua, which further leads to a small A|5| (or small A/?j due to Eq. (l46l) ). As 
expected from the previous discussion, the large |A^J appears in the large 7 region. 
For 7 = 7r there exists a Z2 discrete symmetry in the neutral Higgs potential which 
might cause a dilution between matter and anti-matter bubbles for SCPV-induced 
EWBG (the explicit CP breaking from the fermion sector avoids problems with cos- 
mo logical domain walls [106] ). The vacua degeneracy is lifted as long as 7 deviates 
from 7r. Because the nucleation rates of these two vacua are exponentially biased by 
the difference of their energy densities, a slight shift from ir for 7 can greatly diminish 
the dilution. In this article therefore we will choose 7r/2 < 7 < 71. 

5.3 New Contributions to Electric Dipole Moments 

In this subsection, let us discuss possible new contributions to the EDMs in the SCPV- 
induced baryogenesis (also see [70]). According to the discussions in the previous 
subsection, no CP violation occurs at zero temperature for 7 = 0, it, so there is no 
EDM problem. 

For 7 7^ 0, 7r, the Higgs fields may obtain nonzero phases at zero temperature. 
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However, their contributions to the EDMs are highly suppressed. Among the four 
gauge independent phases, (3% and /3 4 are zero due to the A h and A x terms, so there 
will be no phases entering the h and terms after we do the field redefinitions. The 
only phases are associated with the m 2 SSi m 2 SS2 and rn 2 SiS2 terms. (The phases will 
enter the fermion Yukawa couplings also, but they only lead to overall phases in the 
mass matrices which do not affect the CKM matrix.) Therefore, the CP violation 
will only affect the Higgs sector. 

The Feynman diagrams contributing to the electron and quark EDMs which 
are mediated by the Higgs fields involve two vertices containing the Yukawa couplings, 
and include a mass insertion on the fermion line to flip the chirality. The EDMs will 
therefore be proportional to the cube of Yukawa coupling, making it very small. The 
EDMs are also proportional to the CP violation term which comes from the secluded 
sector and couples very weakly to the ordinary fermions. In particular, if any of the 
three soft mass terms m 2 SSl m 2 SS2 and rn 2 SiS2 vanish, there will be no tree- level CP 
violation. Also, for the parameters that we consider, the imaginary part of tu 2 s s is 
very small. Combining all these effects, new EDM contributions from the Higgs sector 
are very small. Our numerical estimate shows that they are 6 to 7 orders smaller than 
the experimental upper limit. 



6 Wall Thickness and Velocity 

In the early Universe the first order EWPT is realized by nucleating bubbles of the 
broken phase. To calculate the baryon asymmetry, understanding the dynamical 
behavior of these bubbles, such as the nucleation rate, expansion velocity, and bubble 
wall profile, is important but analytically difficult. In this section, we will focus on 
the wall profile and wall velocity. 

6.1 Wall Thickness 

The kinetic Lagrangian for the complex scalar fields can be written as 

i 

where <f>i and 9^ collectively denote the magnitudes and phases of the Higgs fields, 
respectively. For the unitary gauge, we are left with 10 dynamical equations (with 
respect to 6 magnitudes and 4 gauge- invariant phase variables). Once the station- 
ary solutions to these 10 equations are found, we can obtain the relevant physical 
information on the bubble wall, such as the wall thickness and velocity. However, 
these classical field equations are extremely complicated, so we have to solve them 
numerically. 

A useful numerical method was devised in [72, [731 dj- The basic idea [72] is 
that to solve the classical field equations one has to find the field configuration for 
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which 
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-oo c 



(4? 



where E^(z) and E^. (2;) are the relevant field equations under the planar wall approx- 
imation. (Ej3.(z) is rescaled by a factor i so that E^^^z) and E J phase (z) have 
the same dimensions.) This can be achieved by searching for the absolute minimum 
of the action Sa- We need to find an appropriate ansatz satisfying the boundary 
conditions at infinity as a first step. The kink ansatz is especially suitable since it 
satisfies the boundary conditions very well and has a smooth, gradually changing 
behavior as one moves from negative infinity to positive infinity. Therefore, we will 
take the kink ansatz for the independent field variables 

m = W. + W/ + W.-W/ tMh( » )| 

m = (Mi±m + (Mi^Mi tanh( £) 2 , 4 m 

where the index "i" runs over the six complex neutral Higgs fields: H%, H®, S, S±, 
S2, S3, and "t" and "f denote the VEVs in the true and false vacua, respectively. In 
our convention the true vacuum is located at z — > 00 and the false one is located at 
z — > — 00. For simplicity, we assume that there are no off-sets which shift the fields 
against each other. 

However, this assumption does not hold for /3 3 . From the Ah term in Eq. 039p . 
it is easy to see that, as (H^) and (H®) become non-trivial crossing the bubble wall, 
(P 3 ) will be suppressed to zero quickly. So a good ansatz for /3 3 should be 

fc(z)=<^[l-tanh(^)]. (50) 

This is an anti-kink function and a > 5. For simplicity, we set a = 5, so all kink 
ansatzs are mediated by only one parameter - the wall thickness 5. This discussion 
does not contradict the assumption A/3 3 = /3(oo) — /3(— 00) = in Subsection 15.21 
since there we are analyzing the phase change crossing the wall while /5 3 ~ inside 
the wall. 

Figure [5] describes the dependence of the wall thickness on 7 (| < 7 < ir) by 
minimizing the action Sa over S for the parameter sets in Table [3j Again, the three 
curves correspond to the three different m 2 „ values: 2, 4 and 6. These plots tell us 
that the wall thickness becomes larger as 7 increases and reaches its peak value around 
7 = 7r. In the action Sa the kinetic terms and potential slope terms are proportional 
to 5~ 3 and 5, respectively. The wall thickness is a result of the balance between 
them. According to the discussions in Section \5[ large 7 induces large magnitudes 
for A/3j or A#<^, which can increase the weight of the kinetic terms in Sa- The wall 
thickness therefore is an approximately monotonic increasing function with respect to 
7. Numerical simulation shows that the wall thickness varies in the range (3 — 30)T C _1 . 
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Figure 5: Wall thickness 5 vs. 7(7!"). 



6.2 Wall Velocity 

The velocity of the bubble wall also plays a crucial role in EWBG. As a technical 
requirement, the wall velocity cannot be smaller than 0.01. Otherwise, the effect of 
EW sphalerons and Higgs-mediated (in the thin wall regime) decay processes cannot 
be neglected while the CP- violating current is propagating in the false vacuum, and 
the new relaxation effect of EW sphalerons will appear before the wall recaptures the 
left-chiral particles composing the flux. 

The wall velocity of a new-born bubble is mainly mediated by three kinds 
of forces: 1) the driving force (Fj) responsible for the bubbles' expansion; 2) the 
shrinking force (F s ) due to the surface tension of the bubbles (T s ); and 3) the friction 
coming from the plasma (Ff). Though some progress has been reported in the MSSM 
recently [91], a systematic study on the dynamics of the bubble wall is still absent 
in the sMSSM. In the following, we will argue that the wall velocity in these models 
should be comparable to or smaller than that in MSSM or NMSSM under the thin 
wall approximation [92^1 . 

Consider two bubbles from the sMSSM and MSSM, respectively, and assume 
that the energy density differences across the wall, e = V out — V in , are the same for 
both bubbles. In this section, a prime will be used to denote the quantities from the 
sMSSM. The driving force per area acting on the walls is Fd = e and the shrinking 
force is F s = 2T s /r, so we can approximately estimate the critical radii for both kinds 

9 This "thin wall" is different from the one defined in next section in that the comparison scale is 
the bubble radius rather than the mean free path of the scattered particles. 
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of bubbles by employing the equation Fd = F SJ which gives 



r c = 2T s /e, (51) 
r' c = 2T' s /e. (52) 

For T' s = aT s , we have r' c = ar c . Here a ^> 1 since the large VEVs of the Higgs fields 
in the secluded sector will contribute significantly to T' s according to the second term 
of Eq. gTj). 

Now suppose that both bubbles sweep a distance Ar, i.e., the bubble radii 
respectively increase to r c + Ar and r' c + Ar, and that the wall tensions T s and T' s are 
left invariant during this process; then the shrinking forces due to the surface tension 
will become smaller. The net force for both bubbles are changed from zero to 

2T 

AF = F d -F s = e-—^—, (53) 

r c + Ar 

IT' 

AF' = F' d -F> = e--^- (54) 

r' c + Ar 

and 

AF - AF' = 2T S ( — A + — — ) > 0. (55) 

r c + Ar r c + Ar / a 

That is, in the sMSSM the shrinking force decreases more slowly if there is no other 
physics involved (note that the driving forces stay unchanged during this process). 

On the other hand, due to thermal scattering from particles in the hot plasma, 
the bubble wall will suffer friction or thermal pressure Ff during its expansion, where 
Ff is determined by the thermal environment and the wall velocity relative to the 
plasma. Now suppose that bubbles in both models expand at the same velocities, 
microscopically this thermal pressure is mainly induced by momentum transfer of the 
scattered particles, particularly those with 

k-.kf 

m(oo) 3> m(— oo), % ' < 0. (56) 

| fci 1 1 Kf | 

Besides the ones in the MSSM, the sMSSM contains extra light degrees of freedom, 
e.g., the singlets in the Higgs sector and the additional exotic particle species. We 
may well expect that in the sMSSM the thermal pressure is not smaller than that in 
the MSSM given T c = T' c and v w = v' w . Therefore, the wall velocity in the sMSSM 
should be non-relativistic. In this article, we will set v,„ = 0.05. 



7 Non-local Electroweak Baryogenesis 

In EWBG scenarios, the non-local baryogenesis or charge transport mechanism [75] 
refers to the case where particles have CP-violating interactions with a bubble wall. 
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These CP- violating interactions cause an asymmetry in some CP-odd quantum num- 
bers or transport charges other than B carried by particle currents into the false vac- 
uum. In front of the expanding bubble wall, where the thermal equilibrium is broken 
to some extent, the asymmetry of transport charges induces a non- vanishing number 
density for left-chiral particles carrying EW charges. This asymmetry is subsequently 
converted by the B + L- violating EW sphaleron processes, such as %%%t-^ <-> and 
^r^r^r^r 0' m t° an asymmetry in baryon number. Some instant later the wall 
sweeps over the region. The generated net rig is then frozen in the bubble due to 
the exponentially suppressed EW sphaleron effects. There are two typical regimes 
relevant to non-local baryogenesis, depending on the particles under consideration: 
the thin wall regime for 5/ A < 1 and the thick wall region for 5/ A > 1, where 5 is 
the wall thickness and A is the particle's mean free path or propagation length. 

For the thin wall regime [93l IM] , the physical picture is relatively simple since 
we may neglect the damping effect on the CP- violating source caused by the self- 
energy thermal correction. The reflectivity and transmittivity of particle scattering 
can be calculated in quantum mechanics, where the fermions interacting with the 
bubble wall are treated as free particles with space-dependent masses, and the CP- 
violating source can be modeled as a delta function current. 

For the thick wall regime, the effects of additional decay processes in the 
wall cannot be neglected, and hence modeling the CP-violating source will be very 
complicated. Generally, two parallel approaches are available: the semiclassical force 
method [95l [96] and spontaneous baryogenesis [97]. The first one is based on the 
fact that, under the WKB approximation (p ^> S^ 1 ), a CP- violating mass term in 
the Lagrangian can induce a CP- violating shift in gradients of the dispersion relation 
at first order. This CP-violating source then appears in the diffusion equations, 
interpreted as one kind of "semi-classical" source. This approach has been applied to 
the MSSM [90] and the non-minimal NMSSM [15]. However, this approach neglects 
the contributions from the non- WKB region, which may greatly enhance the produced 
baryon asymmetry. The second approach is a little like the one used in the thin wall 
regime: the CP-violating source is the particle current reflected by the bubble wall, 
but the damping and multiple-scattering effects are included. 

With CP-violating sources produced, the particles carrying the CP-odd trans- 
port charges will undergo a diffusion process outside the bubble. Dynamical baryoge- 
nesis requires a slight departure from thermal equilibrium in front of the wall, imply- 
ing that not all particle "decay" processes (i.e., those which affect the charge number 
density) are fast or adiabatic compared to the expanding velocity of the bubble wall. 
If 

DT , . 

2 > 1 (57) 

v 

w 

with T being the decay rate and D the associated diffusion constant, the processes 
will be considered as fast, and vice versa. In our calculations, almost all processes 
included in the diffusion equations are fast, otherwise they will not have enough time 
to significantly affect the charge number density. The only exception is the EW 
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sphaleron process. We assume that they are slow during the diffusion and can be 
neglected in mediating the charge number density. We can therefore deal with the 
diffusion of the particles carrying transport charges and the generation of the net 
baryon number separately. We must then: (1) solve the diffusion equations of the 
particles carrying the transport charges, and find the number density of all left-chiral 
particles til which carry EW gauge charges; (2) insert ni into the baryon diffusion 
equation, to find the generated baryon asymmetry. The description of the particle 
transport in front of the bubble wall by the diffusion equations is well-justified because 
v w «C v s . Here v s = 1/V3 is the sound speed in the plasma. 

In the following, we will study the roles in EWBG played by different particles 
in the sMSSM. We will calculate the contributions from leptons in the thin wall 
regime, and the contributions from squarks, charginos and quarks in the thick wall 
regime by assuming the spontaneous baryogenesis mechanism. 



7.1 The Lepton Contribution — Thin Wall Regime 

In the supersymmetric extension of the SM, leptons generally have a mean free path 
in the thermal plasma much larger than the bubble wall thickness. To consider their 
contribution to EWBG, we can work in the relatively simple thin wall regime [93] .In 
the thin wall regime, the scattering of particles by the bubble wall is described in terms 
of freely propagating particles with space-dependent masses, with the Lagrangian 

C = ^ L i^d^ L + ^ R i^d^ R -m{z)^ L ^ R - m*(z)ip R ip L , (58) 

with the mass term determined by the Yukawa coupling. For charged leptons it is 

m l {z)=y l v l {z)e w ^ . (59) 

Since the Yukawa coupling y T ^> y e ^, we will only study the contributions from r 
leptons. 

The CP- violating current in this case is produced by t r — > tl, where the 
CP-odd transport charge is axial r number 

Q; ifl = Diag(l,-l). (60) 

Explicitly solving the Dirac equation and only considering the contribution from the 
non-WKB momentum subspace (k* < 5" 1 ) as a conservative approximation, one 
obtains the CP-odd reflection asymmetry [531 El] 

m T (oo) 

mJz) 



/oo 
dz cos(2/^)Im 
■oo 



m T (oo)\ 

m r (oo) 2 5A^— — , (61) 

where = tanh</> and tanh(20) = |m T (oo)|/|fcf|. In the second line, the kink ansatz 
defined by Eq. f|49|) and the assumption of A^/2 as perturbative (it turns out that 
| A#i < 0.3 in the parameter region we are interested in) have been used. 
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The axial r current j TL induced by the scattering of leptons is [93], [M] 
r rj 3 k k z 

Jr L = / T^kiUm) - m-k^ARr^ih, k f )Ql L . (62) 

Jkf>m. T V Z7T ) & 

Here ft and // in Eq. (l62l) are the Fermi-Dirac distributions in the wall frame for the 
particles localized in and outside the bubble, respectively, given by 



7 2(! _ Vw K ; 



The chemical potential in the exponential factor is suppressed due to the large value of 
the critical temperature T c . Similarly, due to the exponential suppression for E > T c , 
we may only consider the contribution from the small energy region E < T c and make 
a perturbative expansion with respect to approximately yielding the formula for 
the injected CP-violating current 



-m.Coo^A^ (64) 



where 



is a function of the wall thickness and critical temperature. The first term in Eq. fl65l) 
exactly leads to Eq. (25) in Ref. [91] up to an ansatz-dependent factor after h(S,T c ) 
is integrated out. The second term, which will suppress the injected chiral flux, is 
neglected in Ref. [91]. Note that the CP- violating current j TL is proportional to 
m r (oo) 2 . Given that the associated calculations can be applied to the other two 
families of leptons as well, we are justified to only consider the r. 

While propagating in the false vacuum, for a wall velocity that is not too large 
the behavior of the scattered leptons can be approximately described by diffusion 
processes. The injected current will lead to a "diffusion tail" of the particles in front 
of the moving wall. If this diffusion tail or persistence length £ is much larger than 
the width of the wall, it will be a good approximation to describe the injected current 
as a delta-function-like source. While the current is propagating in the false vacuum, 
leptons will undergo some Higgs-mediated or chiral mixing decay processes which 
may have a negative effect on the CP- violating current amplitude. A reasonable 
assumption is that the characteristic decay time for leptons is much longer than the 
time it takes for the wall to recapture them. We can then neglect lepton decays 
or chiral mixing processes and only consider the left-chiral leptons propagating in 
the false vacuum. Furthermore, we also need to make an assumption that the EW 
sphaleron processes are too slow to alter the density in front of the wall significantly 
and thus can be neglected except for their role in generating baryons. These two 
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assumptions require ^ < 1 for the associated processes. For r leptons, they work 
well if v w > 0.01 [93]. Neglecting these decay processes, the diffusion equation in the 
bubble wall frame can be rewritten as [93] 

D TL n r L ( z ) - v w n' TL (z) + J TL = , (66) 

with 

n TL = dn TL (z)/dz , 

Jr L = -{tJtJ'(z), (67) 

where the persistence length £ TL is given by £ Ti ~ 6D TL (v TL ) [93], with D TL ~ 100 T~ x 
being the related diffusion constant [93] and (v TL ) being the average velocity of the 
injected current relative to the wall; 5'(z) = d5(z)/dz and z is related to the coordi- 
nates in the plasma frame by z = r + v w t. Solving Eq. (IBTj) by combining with the 
boundary conditions 

lim n TL (z) = , (68) 

\z\— >oo 



one obtains 



with 



= j C TL e< iz) z<0 m 
z > ( J 



C TL =6(v TL )j TL . (70) 

The baryon generation is described by the diffusion equation 

D L n" B - v w n' B - n f T ws n L = 0. (71) 

Here Ul is the net number density of left-chiral particles carrying EW charges. Only 
such particles are involved in the EW sphaleron processes and hence able to source 
the baryon generation. Di is effective diffusion constant. In the thin wall region, the 
main contribution to m is from the scattering of r leptons, so we have ul = n TL and 
Dl = D TL . Under the approximation (©(— z) is step function) 

r ws (z) = Q(-z)r ws (-oo) 

= 6Ka 5 w T c e{~z), (72) 
with k, ~ 21 for T c > T EW ~ 100 GeV ESI EH], we find for z > 



n B = n f D ^ C ^ , (73 ) 
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Here the boundary condition ris(— oo) = has been used. Then with s = - ^* c the 
baryon asymmetry is solved to be 



Eq. (I74p tells us that in the thin wall regime of the non-local EWBG scenario, 
three model-dependent physical quantities play critical roles: the variation of the CP 
phase crossing the wall A6\, the wall thickness 5, and its velocity v w . A9\ will bias the 
symmetry between the dynamical behaviors of particles and anti-particles. For the 
perturbative approximation in the non-WKB momentum region, i.e., both A9i and 
momentum are not large (Adi/2 < 1 and fcf < <5 _1 ), the CP-violating current and 
baryon asymmetry are proportional to A9\. It turns out that the contribution from 
the thin wall regime is non-trivial compared to the ones from the thick wall regime. 
The wall configuration is also essential. It measures the size of the region in which 
the CP-violating current is generated, and it provides the p z upper limit or definition 
standard of non-WKB particles. From Eqs. (]64p and (|74|) . the CP- violating current 
and baryon asymmetry are proportional to the wall thickness, and are also affected 
by the wall thickness through the integration upper limit on p z . In addition, l/v w 
measures the time for which the EW sphaleron processes continue before the particles 
in the CP- violating current are recaptured by the bubble wall, so it is natural that 
the baryon produced asymmetry is inversely proportional to v w . 

7.2 The Squark, Chargino and Quark Contributions — Thick 
Wall Regime 

In this subsection we will calculate the baryon asymmetry generated by squarks, 
charginos and quarks by employing the spontaneous baryogenesis mechanism. Com- 
pared to the thin wall regime, there are two new effects which are important for our 
understanding of the interactions between particles and the bubble wall: 

(1) The effect of thermal scattering. This can lead to a self-energy correction 
to the particle propagators and hence suppress the CP- violating source (due to the 
shortened life-time of the particles), and will be taken into account by including the 
imaginary part of the fermion self-energy in the dispersion relation, i.e., 



s 



540 7 > TL )P T£ m T (oo) 2 £Afl 1 M(5,r c )F 
(2tt)V,S*P c 3 



(74) 



[^(k)+*7] 2 = k 2 + m 2 , 
where 7 is the associated damping rate. This immediately gives 



(75) 



3(k)| = i 



(76) 



with 



A = T 



9?(k) 



1 



(77) 



1 
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Since r is the lifetime of the scattered particles, and 3ft(k)/u; is their velocity, A 
defines the propagation length of the scattered particles in the wall. In addition, 
from Eq. (175]) . we have 

u 2 = ^( k ) 2 + m2 + 7 2 (l-^)- (78) 

We will work in the limit j 2 <C !K(k) 2 + m 2 , which gives the relation 

uo 2 = 5K(k) 2 + m 2 . (79) 



(2) The effect of multiple scattering by the wall. This effect can be inter- 
preted as multiple insertions in a perturbation regime. Currently, there are two main 
insertion methods: mass insertion [76] and Higgs insertion [3T] [77]: 

Mass insertion assumes that the tree-level physics is described by kinetic energy 
terms and considers the space-dependent mass terms as a perturbation, which can be 
denoted by 

C = C tree + C pert = K-M, (80) 

where K is kinetic energy terms and M is mass terms. Since tree-level physics in 
the mass insertion regime is a Lorentz-invariant free theory of massless particles, this 
method is effective only if K 3> M. However, due to the Boltzmann factor suppression 
for the injected CP-asymmetric current, the main contribution to the CP-violating 
source comes from the non-WKB momentum region, which usually makes multiple 
mass insertions unstable. 

Compared with mass insertions, Higgs insertion is more extensively assumed 
in the literature. Its starting point is to separate the mass terms into a free part and 
perturbative part, i.e., [31] 

£ = C tree + C pert = (K - M°(z)) - 5M(z, Zj) , (81) 

where z denotes the point where the CP- violating current is calculated, and also 
the mass eigenstates are defined, and Z{ denotes the point where the scattering oc- 
curs. Using these local mass eigenstates as approximate asymptotic states, one can 
calculate the associated S matrix and the CP- violating source. If these local mass 
eigenstates are not the gauge eigenstates in the false vacuum, after the CP-violating 
source is found one needs to go back to the gauge eigenstates by taking an inverse 
unitary transformation to evaluate the diffusion equations in the false vacuum. This 
method is justified for any scattered particle species, but it requires that the Higgs 
field variations be small within a propagation length. This implies A ^ 5. This 
relation adds an upper bound on the particle momentum in the z direction. Numer- 
ically, its effect is weak, since the contributions from the small momentum region to 
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the CP- violating current are more important, due to Boltzmann factor suppression 
in the high momentum region. 

We will use Higgs insertion method. For given momentum and energy, the 
CP- violating current at point "2" is mainly generated in the layer between z and 
z + A. Since A <C 5, at the lowest order we can use a formalism similar to that in the 
thin wall regime to approximately calculate the current. Then the generated baryon 
asymmetry can be solved by embedding this CP- violating source into the diffusion 
equations in the thick wall regime. This approximation is adequate for our purpose, 
i.e., for presenting typical features of EWBG in the sMSSM. 

(I) The squark sector 

The mass square matrix of the up-type squark is 





where the diagonal entries are defined by 




m % + m l + A « l + A ^ > 

™\ + < + A * R + A k > 

yui {hsv x S a W- A hu v 2 e- i9 *) , 
y Ui (hsv 1 e- i ( a+ ^ - A hu v 2 e ie2 ) , 




with 





1 



1 




can be perturbatively expanded at u z" as 
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In the first line, we have made a first order approximation for M~. In the second line, 
"Zj" denotes the space point where the particle scattering occurs. 

According to the bilinear interactions, a CP- violating left-chiral squark current 
jtf is produced through u l R — > u % L , where the CP-odd transport charge is axial top 
number 

Lfl = Diag(l,-l). (86) 

Considering that this is mediated by the Yukawa couplings of the up-type quarks and 
that the Yukawa coupling of the top quark y t = y u . A is much larger than the other 
two, i.e., y Ul ,u 2 , we neglect the contributions from the first two family squarks. After 
a field theory calculation (for details, see Appendix [C]) , we find 

^ R t R ^t L i k i,k f ,z) = — — — — / dz x dz 2 sin((k* -k z f )(zi- z 2 )) 

^hihf v^f J z 

{w 1 (z 1 )w 1 (z 2 )sm(9' 1 (z 1 ) - 9[(z 2 )) 

+w 2 {zi) w 2 (z 2 ) sin(-0 2 (zi) + 9 2 (z 2 )) 

+wi(z 1 )w 2 (z 2 ) sin(9[(zi) + 9 2 (z 2 )) 

+w 2 (z 1 )w 1 (z 2 ) sin(-9 2 (z 1 ) - 9' l (z 2 )) 

-(*! - z) - (z 2 ^ z)}, (87) 



with 



E\ = E) = {hf + ml = (k f f + ml 

E 







= {h? + m 2 iR = 






— _L „* - 


Wi 


(*) 


= y t hs(z)vi(z), 


w 2 


(*) 


= -y t A ht v 2 {z), 




(z) 


= a(z) + 9 1 (z), 




Tt 


= Tt R +Tt L - 



Here the indices "i" and "f" denote the initial and final states of the scattering, 
respectively; and (zj — > z) denotes terms similar to the first four ones in {}, but 
with Zi replaced by z. In this paper, we set the damping rates 7£ = 7£ = 0.1T C ~ 
ct s T c as an approximation, since an analytical calculation for them is still absent. 
ARt R ^ L (ki, kf, z) is sensitive to the A ht soft parameter due to its dependence on w 2 . 
This makes its contribution to EWBG very different for small and large Ah t (EWPT 
is not sensitive to Ah t , since it enters the neutral Higgs effective potential only at loop 
level). Particularly, this contribution can be quadratically enhanced by a large Ah t 
through the w 2 w 2 term. We will discuss more on these issues in the next subsection. 
Using these results, the space-dependent CP- violating current is solved to be 

f d 3 k k z / - \ 

3lM = J kZ>0 J^ys/ f - flU^-W) &Rt R ^ L (ki,k f ,z)Qi. (89) 
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Here fz R (ki) and /z+a( — ^iO are the thermal distributions of in the wall frame at 
"z" and "z + A", respectively. Compared to the one defined in Eq. (1631) . they have a 
factor 3 difference due to QCD color as well as a sign-flip for the second term in the 
denominator. 



(II) The chargino sector 



The chargino mass matrix is 

with Ui = g 2 Vi/y/2 and yt, = hs. It can be perturbatively expanded around z as 

M° ± = M^(z), 

( u 2 (z t )e- ie ^-u 2 (z)e~^\ 

* \ Ul (zi)e- ie ^ - Ul {z)e- i9 ^ fi{ Zi )e ia M - fi{z)e ia ^ J ' { } 

For the parameter values that we are using, |mi,2(z)| <C M 2 ,\/j j (z)\, so we simply 
neglect the off-diagonal elements in M~±. In this sector, the CP- violating Higgsino 

current is produced through W c — > H c , where the CP-odd transport charge is vector 
Higgs number 

Q H - -+--= Diag(0, 1,0,1). (92) 

W+,H+,W ,H toV ' K 1 

After a field theory calculation (for details, see Appendix [C]), we find 

AR Wc ^H C (ki,k f ,z) = —— — — / dzxdz 2 sm((k% - k z f )(z x - z 2 )) 

{fc?fc/M*i)ui(* 2 ) sm(9i(zi) - 9i(z 2 )) 
+u 2 (z 1 )u 2 (z 2 ) sin(-d 2 (zi) + 9 2 (z 2 ))} + 
M 2 n(z)[u 1 (z 1 )u2(z 2 ) sin(0i(zi) + 9 2 {z 2 ) + a(z)) 
+u 2 (z 1 )u 1 (z 2 ) sin(-9 2 (z 1 ) - 9i(z 2 ) - a(z))] 
-( Zl ^ z )-(z 2 ^ z )}, (93) 

where 

E} = E) = {hf + Ml = (k f ) 2 + fi(z)\ 
v* - *L v * - */ 

7x c = 7w<= + 7h=- (94) 



10 In the chargino sector, the CP-odd transport charge can also be axial Higgs number. But this 
source usually is suppressed by a fast Higgsino violating process corresponding to the h(S)hih 2 
term in the Lagrangian |31j . We will therefore not consider it here. 
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For the damping rates, we will use the results in the MSSM, which are given by 
Ifjc = 0.025T C and 7^ = 0.065T C [lOOj . The CP-violating current is solved to be 



d% k) 
kf>o W?E~f 

(jf\K) - /2 C A(, xC )(-^)) AR Wc ^ c (k t , k f , z)Ql +x , (95) 

where f^°(kf) and /]^a( — M) are the thermal distributions of the charged gaugino 
W c in the wall frame, which have a factor 2 difference from that defined in Eq. (1631) 
due to W c = Wt + W C R . 

(Ill) The quark sector 

The top quark mass is 

mt = h t v 2 e ld2 . (96) 
which can be perturbatively expanded around z as 

m\ = m t (z), 

5m t = h t v 2 ( Zl )e ie ^ - h t v 2 (z)e id ^ z \ (97) 

In this sector, the CP- violating current is produced by t R — > t^, where the CP-odd 
transport charge is axial top number 

<4 K = Diag(l,-l). (98) 

Then (for details, see Appendix [C]) 

AR tR ^ tL (k h k f ,z) = — — — / dz^sin^k* - k z f )(z l - z 2 )) 

HiiEjf V i V f Jz 

{h 2 t k z k z v 2 { Zl )v 2 {z 2 ) sm(9 2 ( Zl ) - 6 2 {z 2 )) 
-( Zl ^ z )-(z 2 ^ z )}, (99) 



where 



Et = E) = {KY + K°(z)|' = (k f Y + \m° t (z)[ 

k z 



It = lt R + lt L - (100) 

Again, we approximate the damping rates by % L = jt R — 0.1T C ~ a s T c . The trace 
in the S matrices of the chargino and top quark scattering gives an extra factor "4" , 



37 



so the factor "4" in the denominator of Eq. (jC.llj) disappears in both Eq. fl93|) and 
Eq. (199)) . The CP- violating current is solved to be 



3t L \ z ) 



d 3 ki k 



f 



f 



f*(kX) - . * fiUtU-K) &Rt R ^t L (ki, fe,^)Q* /ioi) 



fcf 2 - Am 2 

where fz R (k*) and /z+a( — ^f) are the thermal distributions of the right chiral top 
quark in the wall frame, which have a factor 3 difference from that defined in Eq. (163)) 
due to color, 

Am 2 = m 2 + A(7t)) -m 2 t {z) (102) 
is the difference of the quark mass square at 2 and 2 + A (7$), and 



A^ = ^m2(oo)-7n?(z) (103) 

is the momentum cutoff. The factor in front of ft L (—kf) in Eq. (11011) is due to the 
large mass difference of the scattered particles inside and outside the bubble. For r 
leptons, stops and charginos, this difference is small compared to T c , so this factor is 
approximately equal to 1. 



Before moving to the calculation of the produced baryon asymmetry, we have 
some comments on the CP-violating currents obtained from the lowest Higgs inser- 
tion. First, according to the discussions in Subsection I5.2[ the single Higgs phases 9^ 
are not physical since their values depend on two decoupled degrees of freedom: A 
and B. Indeed, though working in the 9^ basis, from Eq. ( 1ST)) . Eq. ( 193)) and Eq. ( 1991) 
we see that all of these currents depend on some special linear combinations of 9^ , 
where their dependences on A and B are cancelled and they are only sensitive to 
(3i(z). Second, for the CP- violating currents from the stop and the chargino sectors, 
the terms oc uj\uji and oc u\Ui will be reduced to the v\dv-i — v-idv\ type sources in the 
MSSM if 9\fi{z) = and a(z) = constant. Though in this limit these contributions 
are trivial for dtan/3 = 0, they are not in Eq. (1ST)) and Eq. (J93J) due to the space- 
dependence of #1,2 in the wall. Third, the terms oc ViVi in Eq. (1ST)) and Eq. (1931) are 
zero in the limit 9i^{z) = and a(z) = constant. They therefore never appear in 
the MSSM currents and denote a new-type of CP-violating current in the stop and 
chargino sector. This is also the only type of CP-violating current at the leading 
order in the quark sector. The vector Higgsino current of the Vidv 2 + v 2 dvi type 
does not appear in our calculations since it denotes a resummation effect from higher 
order Higgs insertions [TTJ [3H [32]. Finally, we emphasize that all of these features 
are results of SCPV-driven EWBG, so they are not sensitive to the concrete embed- 
dings of the sMSSM. (For different embeddings, the U(l)y and U(l)' gauge charges 
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in Eq. (]43p and hence the formulas in Eq. (I4ip are different. However, we did not use 
these results when deriving the formulas of the CP- violating currents.) 

The field theory calculations used in the thick wall regime cannot be directly 
applied to the thin wall regime. In the thin wall regime, the main contribution is from 
an interference effect between a CP and a non-CP term (see the Appendix of [94|). 
The non-CP term is induced by the effect that the fundamental solutions inside 
and outside the bubble are different due to the mixing between left- and right-chiral 
fermions. In the quantum mechanics calculations of the reflection and transmission 
coefficients, this effect is explicitly included. Its strength can be measured by the 
parameter which is defined in Subsection 17.11 In the thick wall regime, this effect 
is neglected in the field theory calculations. Instead we assume the same asymptotic 
states between the two sides of a layer inside the wall, which greatly simplifies the cal- 
culations. This assumption is justified by the fact that t$ is approximately estimated 
to be 

Am(z) A 

U * ~ 1 M ~ « 1 (104) 
Am 

in the thick wall limit. Here Am(z) = m{z + A) — m(z) is the mass difference of the 
incident particles between the two sides of the layer, Am = m(oo) — m(— oo) is their 
mass variation crossing the wall, and A is the thickness of the layer (see Eq. (177)) ). 
The main contribution in the thick wall regime therefore is from an interference effect 
of the CP- violating terms. Actually, if we apply the field theory calculation to leptons 
(similar to what we did for quarks), we will get their subleading contribution, which 
is given in the Appendix of [94"] . 



While these injected CP-violating Higgsino or squark currents propagate in 
the false vacuum, the original thermal equilibrium among different particle species 
is broken and these species may obtain nonvanishing net number densities due to 
the associated decay processes. In the thick wall regime, the main decay processes 
include [75] : (1) the Yukawa interaction T yT corresponding to all supersymmetric and 
soft trilinear interactions arising from the UtH2QlTr term in the superpotentiaj"! 
(Ql is the left-chiral superfield of the third family quarks and Tr is the right-chiral 
superfield of the top quark); (2) the Higgsino violating process T M corresponding to 
the h(S)h\h2 term in the Lagrangian and supergauge interaction T g ; (3) the axial top 
number violating process T m and the Higgs- violating process T h corresponding to top 
quark mass effects and Higgs self interactions, respectively; (4) the strong sphaleron 
process T ss [1U1] and weak sphaleron process F ws . T m and T h are suppressed in the 
false vacuum and T ws is suppressed in the true vacuum (for a strong enough first 
order EWPT); all of the others can occur both inside and outside of the bubble. For 
simplicity, we neglect all processes with a singlet component involved in the current 
propagation. In the concrete calculations, further simplifications are also made |76j . 



11 As a first order approximation, we neglect the tau and bottom Yukawa interactions. For relevant 
discussion, see [102] . 
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(a) We assume fast enough T g which thermally equilibrates n q and riq] thus we can 
describe the system by the densities uq l = n qL +riq L , n TR = n tR + n^ R , n Hl — n^+ng 
and n H2 = n h2 + rijj . (b) We assume fast enough T^; then we have n Hl = n H2 and 
hence are able to simplify the two quantities and n# 2 into one tir = nn 1 + nn 2 ■ 
(c) We assume slow T ws , which allows us to ignore leptons for the particle diffusion 
and the til generation. Under these approximations, the diffusion equations in the 
wall frame ar^l [76] 



with 
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(105) 



(106) 
(107) 



(108) 



where (vx t fj c ) is average velocity of the scattered particles and {kQ L , kT R , kB R , kn} 
are statistical factors. 

For light particles compared to T c , their statistical factor is 2 for bosons and 
1 for fermions. If these particles are very heavy, then their statistical factors are 
suppressed exponentially. Explicitly, this suppression factor is (e.g., see [103j ) 



12 



m \ 3 / 2 



2nT c 



exp 



-m 



Tr 



(109) 



with m being the particle mass. In this paper, we assume that all charginos and 



A = m% 

while keeping the other squark soft mass squares to be 25 (see Table [3]). This soft 



neutralinos are light for simplicity. As for squark masses, we take m - = mt 3 

V3 U R 



12 These equations arc only valid in the limit of non-relativistic wall velocity. This is also true for 
Eq. (671). 
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mass pattern is phenomenologically favored due to the large beta functions for mj~ 
and m? 3 . According to Eq. ( 11091) . we immediately fincf^l 



k QL ^10, k TR ^5, k BR ^3, k H &12, (110) 



and 



-5 23 
D = —D 0j + — D H , 
28 Ql 28 ' 

r = — (r m + r fe ), 

336 

23 

J = ^(\ + Jt L + J^)- (HI) 

Here D, T, and J are the effective diffusion constant, decay rate and CP-violating 
current (for their definition, see [76]). It has been implicitly assumed that the other 
Higgs singlets are not thermally equilibrated with the EW charged ones. Using Dh ~ 
110/T C and Dq l ~ 6/T c [93], one finds D ~ 91/T C . This result further constrains the 
wall velocity: slow weak sphaleron processes require 

DT ws /v 2 w < 1, (112) 

or 

v w > 0.02. (113) 

As for f , using the estimation for T m (z) + Th(z) with tan/3 > 1 in [76], we have 

f (oo) « 10~ 2 T C . (114) 

To solve Eq. (U05H107j) . we further assume that T ss and T yr are large, and 
approximate 

e(z)J(z) -> J(z), 

e(2)r(oo) -> r(*), (115) 

following [3TJ[76]. Then using the boundary conditions ^(±00) = 0, we find 

^=^+^1 = ^/^ (116) 



for z < 0, where 



& 5 1 



L 



46 DA_ 



/■oo 

/ rfC^(C)e" CA+ , (117) 
^0 



1:5 



As an approximation, we use their soft masses as the squarks' physical masses. 
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vs. 7(7r) for cases a, b, and c denned in Table [3l 



with 



4T(oo)D } « 1(T 2 T C . 



118) 



This solution is the same as the MSSM one obtained in [211 EH] except the coefficient, 
which is caused by the difference of the {kq L) kT R , kB R ,kn} values. Solving the diffu- 
sion equation for rig, which has exactly the same form as Eq. (ITT]) with Dl — > D and 
Cl — » Cq l , we obtain the baryon number density for the z > region 



n B 



-rifDT ws CQ L 



;ii9) 



and finally 



n B 



45 15T,, 



27r 2 ^T c 3 46^A + y 



dC J(C) e 



7.3 Numerical Results of Non-local EWBG 



(120) 



According to the previous sections, the explicit CP phase, 7, plays a crucial role in 
mediating the EWBG. Both the spontaneous CP phases and the bubble wall physics 
show a strong dependence on it. There should be a similar 7 dependence, therefore, 
for the generated baryon asymmetry. 

Let us first consider the 7 dependence of v(T c )/T c , the measure of the EWPT 
strength. Some related issues on EWBG have been discussed in Section |4] where no 
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Figure 7: — vs. 7(71") for v w = 0.05, in case a. There is a sign flip for ns/s in 
the region denoted by black-dashed line for the top quark curve. 




Figure 8: ^f- vs. 7(71") for v w = 0.05, in case b. There is a sign flip for ub/s in 
the region denoted by black-dashed line for the top quark curve. The 
stop contribution is smaller than the scale of this figure. 
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Figure 9: — vs. 7(71") for v w = 0.05, in case c. There is a sign flip for ub/s in 
the region denoted by black-dashed line for the top quark curve. The 
stop contribution is smaller than the scale of this figure. 

ECPV is implicitly assumed. As expected, v(T c )/T c shows a strong dependence on 
the soft SUSY breaking parameters m 2 - , m 2 3 and A h . In Figure [6], we illustrate the 

7 dependence of v(T c )/T c . The weak 7 dependence is because 7 is from the secluded 
sector of the neutral Higgs potential, which has no direct coupling to H® and H®. 
For all of the three b and c, the first order EWPT is strong enough for 

implementing the EWBG. 

In Figure [7}l9l we show the 7 dependence of the baryon asymmetry \ub/s\ 
produced by r leptons, top squarks, charginos and top quarkp. We set v w = 0.05. 
The produced baryon asymmetry is sensitive to 7. This is easy to understand since 
both the magnitude of SCPV and the bubble physics are sensitive to 7. Within 
theoretical uncertainties, the observed value %/s = (8.82 ± 0.23) x 10 -11 can be 
explained for a large range of 7 values. Top quarks play a significant role in all 
three cases because of the top Yukawa coupling enhancement. In addition, due to the 
large change of the top quark mass crossing the bubble wall, the cancellation between 
fl R (k*) and /f+ A (— kf) in Eq. (11011) is relatively weak. The stop contribution depends 
on the magnitude of the soft parameter Ah t . For a small Ah t , it is dominated by the 
W1W1 ~ s 2 v\ term in Eq. flH7|) . As m 2 increases, v\ and s usually becomes smaller 

d 

14 Wc did not consider the contribution from neutralinos. Since the relevant CP-violating sources 
and the interactions with the bubble wall are similar for neutralinos and charginos, their contri- 
butions are expected to be comparable with those from charginos. For relevant discussions in the 
MSSM, e.g., see 
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and smaller. The stop contribution is then quartically suppressed. This explains 
why the stop contributions are smaller in case b and c than in case a, and also 
explains why the stop contribution is small for a small m x o in Figure [15] of the next 
section (for all figures in this article, we chose Ah t = 0.1). For a large Ah t , the stop 
contribution is dominated by the W2W2 ~ A\v\ term in Eq. ( IBTj) . It is not sensitive to 
the soft parameter m 2 rr and can be comparable with or even larger than the chargino 
contribution. Numerical results show that the stop contribution becomes comparable 
with the chargino one as A^ t increases to ~ T c and with the top quark one as A^ t 
increases to ~ 4T C . As pointed out in the previous section, the f 2 f 2 -type CP-violating 
current is absent at the leading order in the MSSM. Both the large contribution by top 
quarks and the enhancement effect of the stop contribution by a large Ah t therefore 
are absent in the MSSM. 

In addition to 7, there are two other important factors which can influence 
the produced baryon asymmetry. The first is v w . One can see this from Eq. (1741) 
and Eq. (I119p . A larger baryon asymmetry can be generated if the bubble wall is 
expanding at a more non-relativistic velocity; and vice versa. The second is the stop 
masses. In this paper we assume the stop soft mass I^qJ = \ m u R \ to be heavier 
than T c , even though they are lighter than those of the other squarks and sleptons. 
If they (or more precisely, the stop physical masses) are light compared to T c , then 
the associated statistical factors will be modified to [76] 

k QL ^18, k TR ^9, k BR ^3, fc H «12, (121) 

which gives an additional factor 3 in the numerator of Eq. (I120|) . 

8 Neutralino CDM and Non-local EWBG 

Besides the baryon asymmetry, another mystery is dark matter (DM), which has been 
measured precisely by astrophysical and cosmological probes |lU4j : 

tt DM h 2 = 0.106 ±0.008. (122) 

Given that supersymmetry with i?-parity conservation could have a neutral LSP as 
a good candidate for the latter, e.g., the lightest neutralino x\i one naturally wants 
to know whether or not these two cosmological problems can be solved in the same 
framework. 

To answer this question, we will not scan the whole space, since it is large. 
Instead, we vary parameters to which the neutralino DM relic density or m x o is 
sensitive. From the discussions in the previous sections, m 2 can be such a parameter. 
Explicitly, we use parameter values given in Tabled except that we vary m 2 between 
10 4 ~ 5 x 10 5 GeV 2 and set 7 = 0.957T. The relation between m v o and m 2 is 
illustrated in Figure [TU] As m 2 „ increases, m v o decreases from more than 90 GeV to 
35 GeV. This is because m x o ~ hvi in the limit of large (5*1^,3) and v\ <ti %h ~ (S) (see 
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lower bounds on the DM relic density measured by astrophysical and 
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Figure 14: 5 vs. mo. 
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Eqs. f lA~25llA~27l) ) while a lar ger Tfijj implies a smaller v\. This feature is important 
for the DM relic density, since the Z pole region is covered. As is well-known, the relic 
density of neutralino CDM is usually overproduced. One way to reduce the abundance 
is from the Z-pole region, where we have rrtz ~ 2m x o and the annihilation of two 
Xi through a Z boson is greatly enhanced. In this region, therefore, we can easily 
obtain an acceptable value for the Xi renc density. As a first-order approximation, 
we calculate the Xi renc density today, only counting in the processes mediated by 
gauge bosons (for general discussions on the relevant calculations, e.g., see |105j ). 
The results are shown in Figure [HJ 

Let us come back to EWBG for this parameter region. With m x o (or equiv- 
alently rrvL ) as a varying parameter, the strength of the EWPT, the magnitude of 

d 

CP-violation, the thickness of the bubble wall and the produced baryon asymmetry 
by different particle species are shown in Figures [T2lfT5l From Figure [HJ it is easy to 
see that the measure of the EWPT strength v(T c )/T c is always larger than 1.1. As 
m y o increases or equivalently m 2 „ decreases, its value increases. This can be under- 

stood because a smaller scaled mL (the real m 2 H o value is also smaller) usually leads 

d d 

to a larger scaled V\(T C ) in the true vacuum, while keeping the scaled V2(T C ) and T c 
almost invariant (T c ~ 1.4 in our simulation). In Figure [T3l we see that A^i^ and 
Act are close to zero around m x o « 40 GeV, which explains the small wall thickness 
in Figure [TH (recall that a large phase variation crossing the wall usually leads to a 
small wall thickness (see discussion in Subsection 16. ip ) and a sign flip of ub/s for 
the same m x o value in Figure [HI Another sign flip for the top curve in Figure Q3] 
is related to the ones shown in Figures EHH1 which is a result of balance between the 
contributions from the particle reflection and transmission. Finally, for most region 
of m x o ~ 35 — 90 GeV a strong enough first order EWPT and a reasonable baryon 
asymmetry can be obtained, and in the Z pole region the right Xi re li c density can be 
produced (see Figure [T5l) . The two cosmological mysteries, baryon asymmetry and 
DM, can therefore be explained simultaneously. 

9 Some Other Cosmological Effects 
9.1 Superconducting Cosmic Strings 

There may exist phase transitions in the early Universe due to gauge symmetry break- 
ing, and then topological defects, such as monopoles and/or cosmic strings, etc, can 
be produced |106l 1107] . If the gauge group G is broken down to a subgroup H, we 
will have monopoles and cosmic strings if the homotopy group 7T2{G/H) ^ I and 
7ii(G/H) = Z 7^ /, respectively, where Z is the group of integers. In particular, 
the cosmic strings may be superconducting because of the presence of a charged field 
condensate in the core of the string |108j . The gauge symmetry in the simplest cosmic 
string model is U(l) x U(l)' where the U(l)' gauge symmetry is broken via the Higgs 
mechanism |108j . Note that 7Ti(?7(l) / ) = Z, so we can have a topologically stable 
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cosmic string solution. Although the U(l) symmetry is respected by the vacuum, it 
can be broken in the string core and then give rise to a bosonic condensate. Moreover, 
there may exist fermionic superconducting strings that have fermionic zero modes in 
its core [108] . These fermions can have non-zero U(l) charges, and arise from Yukawa 
couplings to the scalar fields. In the non-supersymmetric models, we only need to 
introduce two complex scalar fields: one is charged under U(l) while the other is 
charged under U(l)' |108j . However, in order to have the bosonic superconductivity 
for cosmic strings in the supersymmetric U(l) x U(l)' model, one needs to intro- 
duce at least five chiral superfields |109j . Interestingly, one can realize the fermionic 
superconductivity in the simplest supersymmetric Abelian Higgs models [110] . 

In the sMSSM, we have the fermionic superconductivity, and supersymmetry 
is broken in the core of the cosmic strings (supersymmetry is not broken outside 
of the string core) which is similar to that in Ref. [110] . However, we do not have 
the bosonic superconductivity since all the Higgs fields are charged under the U(l)' 
gauge symmetry. Because the U(l)' gauge symmetry is broken around the TeV scale, 
the string tension T§ should be around the TeV scale as well. The gravitational 
effect of cosmic strings with T|/Mpj ~ 10 -30 is extremely weak and is not within 
foreseeable detection capabilities. Interestingly, we can search for the signatures of 
particle emissions, such as positrons or other decaying particles produced where the 
current on the string quenches. Similar to the discussions in Ref. [Ill] , which also 
has strings around 1 TeV, the sMSSM might explain the observed 511 keV emission 
from the electron-positron annihilation in the Galactic bulge |112l 11131 1114j . 

9.2 Gravitational Waves 

There are a few experiments under consideration that might detect for the first time a 
stochastic background of GWs, for example, the planned Laser Interferometer Space 
Antenna (LISA) [115] , and the Big Bang Observer (BBO) [116j . which is the follow- 
on mission to LISA. One of the well-motivated and convincing sources for GWs is 
inflation, and the discovery of a relic gravitational background could be a smoking 
gun signal from inflation. Another source of the stochastic GWs are strong first-order 
phase transitions. In contrast with the inflationary spectrum, the spectrum from a 
phase transition is not flat and has a characteristic peak related to the temperature 
at which the phase transition took place. Especially, LISA and BBO will have fairly 
good sensitivity at the frequencies that coincide with the redshifted spectrum of GWs 
produced during an electroweak phase transition at temperatures around 100 GeV. 
Thus, the GW observations will give us information about EWBG since the strong 
first-order electroweak phase transition is a necessary requirement. 

There are two distinct mechanisms that can generate the GWs during a first- 
order electroweak phase transition: the bubble collisions jl 171 11181 11191 H^Oj and 
the turbulent motion of the primordial plasma [Ml ESJ EM The GWs 

from the electroweak phase transition have been studied in a model independent way 
in Refs. |126j and |127j . The results can be presented as a function of two main 
parameters: the typical size of the colliding bubbles and the available energy. It 
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was concluded that a sufficiently strong electroweak phase transition could lead to an 
observable GW signal at LISA. However, for a definite model, the bubble size and the 
available energy are correlated. Thus, it is still very important to study the concrete 
models in details. 

It was found previously that in the MSSM the produced amount of GWs is 
orders of magnitude below the LISA sensitivity. The situation in the NMSSM seems 
much more promising because the trilinear term SHdH u can naturally lead to a much 
stronger electroweak phase transition |128j . Recently, it was pointed out that in 
Ref. [128j a rough approximation was used to calculate the bubble configurations, 
which overestimated the strength of the phase transition and the GW signal [129j . 
Moreover, in Refs. [129[ 1130] the nMSSM, which solves the /x problem and domain 
wall problem simultaneously, was considered. It was shown that the GW signals 
can be detected at BBO for part of parameter space, but cannot be observed at 
LISA pQlfTSO] . 

In the sMSSM the strong first-order electroweak phase transitions arises from 
the trilinear term SH<iH u in the superpotential, which is similar to the NMSSM 
and nMSSM. The GW signals might therefore be detectable at BBO for part of the 
parameter space, but cannot be observed at LISA [129[ 1130] . As a side comment, 
the phase transition for the U(l)' gauge symmetry breaking may not be strong first 
order since the coupling A in the trilinear term XS1S2S3 is small, and one might not 
be able to observe the GWs from the U(l)' gauge symmetry breaking. 

10 Discussion and Conclusions 

In this paper we discussed EWBG and its correlation with the neutralino CDM in the 
sMSSM. We first constructed two anomaly-free models. In Model I, we embedded the 
SU (3)c x 577(2) l*.U (l)y x U(l)' gauge symmetry into a larger gauge symmetry, E 6 , 
whose representations are anomaly-free. We considered three families of SM fermions 
and one pair of Higgs doublets arising from three E$ fundamental representation 
27s. To include the SM singlet Higgs particles, achieve gauge coupling unification 
and cancel the U(l)' anomaly, we assumed that three pairs of vector-like SM singlets 
and one pair of vector-like Higgs doublets from two pairs of 27 and 27* are light 
while the other particles in the two pairs of 27 and 27* are absent or very heavy. In 
Model II, we calculated the anomaly-free conditions, and added the minimal number 
of exotic particles to cancel the U(l)' anomalies. We required that all the particles 
have rational U{1)' charges. Although we had minimal exotic particles, we lost the 
gauge coupling unification. We also presented the general superpotential and the 
supersymmetry breaking soft terms in both models. 

We discussed the one-loop effective potential at finite temperature in the 
't Hooft-Landau gauge and in the MS-scheme. We showed that there exists a strong 
enough first order EWPT because of the large trilinear term AhhSH u Hd in the tree- 
level Higgs potential. Unlike the MSSM, the stop masses can be very heavy. The 
EWPT features in both models are quite similar because the exotic particles' contri- 
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butions to the one-loop effective potential at finite temperature are suppressed due 
to their heavy masses. We therefore restricted the detailed analysis to Model I. 

In the early Universe the first order EWPT is realized by nucleating bubbles of 
the broken phase. Crossing the bubble walls, the VEVs of the Higgs fields (including 
their magnitudes and phases) are space-dependent. Because the dynamical behaviors 
of these bubbles, such as wall profile and expansion velocity, can have important 
influences on the production of the baryon asymmetry, we also studied their physics in 
detail. Numerical results show that the bubble wall thickness varies from 3 ~ 30 T~ x 
as a monotonically increasing function of the phase changes of the Higgs fields. We 
also argued that the wall velocity in the sMSSM cannot be larger than that in the 
MSSM under the same phase transition condition and thus should be non-relativistic. 
This fact implies that EW sphaleron processes have more time to occur and hence 
will enhance the final baryon asymmetry. 

We also discussed possible CP violation introduced by the extended Higgs 
sector. Unlike the MSSM, where there is no CP-violation at tree level in the Higgs 
sector, the CP symmetry can be broken spontaneously as well as explicitly. With 
loop corrections included, there may coexist vacua at finite temperature with broken 
and unbroken EW symmetry. The values of these spontaneous CP phases usually are 
different in these vacua. In our work, the SCPV provides a direct source for baryo- 
genesis while its magnitude is mediated by an explicit one from the secluded sector. 
These new CP sources do not introduce significant new contributions to EDMs. Af- 
ter proper field redefinitions, the CP- violation phases will only appear in the Higgs 
mass matrix through soft masses associated with singlet components. Numerical re- 
sults show that, for typical parameter values, their contributions to EDMs will be 
about six or seven orders smaller than the experimental upper limits. These contri- 
butions disappear completely in the limit with a trivial explicit CP phase, where the 
spontaneous CP phases are trivial in the true vacuum but not where EWBG occurs. 

Subsequently, we systematically studied non-local EWBG in both the thin wall 
and thick wall regimes. We calculated the contribution from leptons in the thin wall 
regime, and the ones from squarks, charginos and quarks in the thick wall regime in 
terms of the profile of bubble wall. For leptons and quarks, the induced CP-violating 
currents outside the bubbles are proportional to Yukawa couplings squared, so we 
only considered the contributions from r leptons and top quarks. We found that the 
CP-violating currents for stops and charginos are very different from those obtained 
in the MSSM. Due to the space-dependence of the relevant CP phases, they do not 
require a variation of tan f3 in the bubble wall to have a non-trivial structure at the 
lowest order of Higgs insertion. In addition, there exists a new type of CP- violating 
current at the leading order which can contribute to the generation of the baryon 
asymmetry. The CP-violating current of the new type has important influence on 
the EWBG. First, the stop contribution can be quadratically enhanced by a large A^ t . 
Second, in addition to r leptons, top squarks and charginos, top quarks can also play 
a significant role in the EWBG. We emphasize that all of these features are results 
of the SCPV-driven EWBG, so they are not sensitive to the concrete embeddings of 
the sMSSM. Numerical results show that the produced baryon asymmetry is large 
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enough to explain the cosmological observations. 

After that, we studied the correlation between EWBG and neutralino CDM in 
the sMSSM. Though we did not scan the whole parameter space, we found that there 
exists a region where a strong enough first order EWPT, large CP phase variations 
crossing the bubble wall, a reasonable baryon asymmetry, as well as an acceptable 
neutralino LSP relic density can be achieved simultaneously. We also commented on 
possible cosmological signals of the model: superconducting cosmic strings and GWs. 
Particle emission from the decays of cosmic strings and GWs from EWPT could be 
observed within the foreseeable future. 

The secluded U(l) '-breaking sector plays an essential role. It induces ECPV 
and SCPV at tree level in the Higgs sector, helping avoid significant new contributions 
to the EDMs which may hinder a successful EWBG in supersymmetric models, and 
also influences the bubble wall physics. It provides more degrees of freedom to drive 
the realization of EWBG, and makes the EWBG physics in the sMSSM very different 
from those in the MSSM and NMSSM. These differences make it interesting to further 
study the relevant collider implications for the LHC [01], which we will leave to future 
exploration. 
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A Mass Matrices for the Particles in Model I 



For the Higgs mass matrices, here we only present the results in the limit m| s = 0, 
where there is no CP violation. It is not difficult to extend them to the general case 
with non-zero m 2 s s and CP violation, as considered in the main text. Let us define 

«i = <flS>, v 2 =(H° u ), tan/5=^, (A.l) 

s=(S) , Si = (Si) . (A.2) 

Then the mass-square matrix for the CP-odd neutral Higgs particles in the basis 
[H? = Im(H° d ), H°J, S°\ S™, S?, S*} is 

K>=(°r£ C c A °) , (A.3) 



where 



and 



C A o S A ° 



P 2 Hd A h hs A h hv 2 

A o = | A h hs p 2 Hu A h h Vl | , (A.4) 

A h hv 2 A h hv x Pi 

p 2 Si A x Xs 3 A x Xs 2 ' 

S A o = ( A x \s 3 Pi A\\si | , (A.5) 

A x \s 2 A x \si p 2 3 

0' 

C A o = | | , (A.6) 

-m 2 SSi -m 2 SS2 



Pk = m 2 Hd + h 2 (v 2 + s 2 ) + ^(v 2 - v 2 2 ) + g 2 ,Q Hd A , (A.7) 



Pk = m 2 Hu + h 2 {v{ + s 2 ) + ^{v\ - v 2 ) + g 2 z ,Q Hu A , (A.8) 



P 2 s = m 2 s + h 2 (v 2 + v\) + g 2 z ,Q s A , (A.9) 
Pi = m% + A 2 s 2 + g 2 z ,QsA , (A.10) 

3& 
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where 

3 

A = Q s s 2 + Q Hd v\ + Q Hu v 2 2 + Y^QsA ■ (A.ll) 



8=1 



Similarly, in the basis {Hf = Re(H%), #° r , S 0r , S?, S 2 0r , 3^}, the mass-square 
matrix for the CP-even neutral Higgs particles is 



where 



4 d K H d ,H u ^H d ,S 

Oho = ( K Hd , Hu 4„ k Hu)S I , (A. 13) 

K>H d ,S K H U ,S 4 



4i K S!,S 2 K s u s 3 
Sh° = I «Si,S a 4 2 K S a ,S 3 | , ( A - 14 ) 
^Si,S 3 K S 2 ,S 3 4 3 



C H o = | Kf^.Sx K^ ui 5 2 K Hu ,5 3 | , (A. 15) 

. ksa + m 2 SSl n s ,s 2 + m\s 2 ^s,s 3 



and 



4 d = 2 + ^'^-) ^ + m n* + ^ + §2) + " ^ + & ,Qh ' A ' (A ' 16) 

4„ = 2 (t + ^ + <. + h \< + * 2 ) + t ( ^ 2 - * + <& Q *« A ' (A - 17) 

4 = 2^|,Q|s 2 + m% + /i> 2 + *; 2 ) + <7§,Q S A , (A.18) 



4 = Igl'QlA + m| + A 2 s) + glQsA , (A.19) 



G 2 

/^ d)jF j M = 2 ( /i 2 - — + g 2 z ,Q Hd Q Hu ) v lV2 - A h hs , 
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refits = 2(h 2 + g 2 z ,Q Hl Qs) vis - \eij\A h hvj , 

KHi,Sj = Zgl'QmQsjViSj , n s , Sl = 2g 2 z ,QsQ Sl sSi , 

Ks^Sj = 2(A 2 + gl'QsiQs^iSj - \e ijk \A x Xs k . 
The charged Higgs mass is 



A/f 2 - M 2 -4- ^ Ah ^ s 
Mh± ~ Mw + ^2P 



-h\vl + vl) 



(A.21) 
(A.22) 
(A.23) 

(A.24) 



where = g\{v\ + vf)/2. 

In the basis {B',B, W®, H®, H®, S, Si, S 2 , S 3 }, the neutralino mass matrix is 



M^o(6,6) M^o(6,3) 
M^o(6,3) T M^o(3,3) 



(A.25) 



where 



M^o(6,6) 



/M> 








pi* 


pi* pi* 





Mi 







73^°* 








M 2 




-7I^ * 













pi* 

1 HP 








hH° d 










hH° u 


hH° d 



(A.26) 



and 



XS 3 XS 2 
M^o(3,3) = | XS 3 A5i | 
XS 2 XSt 



(A.27) 



where = \[2gz>Q<t>§\ and M{, Mi and M 2 are the gaugino masses for 17(1)', £7(1) 
and 5C7 (2) L , respectively. The first row of M^o (6, 3) is given by (r^ T* S2 P^J , while 
the other entries are zero. 

The chargino mass matrix is 



M, ± 



M 2 f 2 H° u * 
f 2 H° d * hS 



(A.28) 



In contrast to the Higgs mass matrix, the neutralino and chargino mass matrices as 
well as the other ones in the following, are written in a general background with CP 
violation. 
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In the basis Wj, Z^), the mass matrix for the neutral gauge bosons is 

M B,W 3 ,Z' M B,W 3 ,Z< M B 3 W 3 ,Z> 

M B \w 3 ,z' M B 2 W 3 !Z , M^ w3zl 



where 



M B \ w3 , z , = \gl[\Hl\*+\H^\, (A.30) 

M B 2 W3 , Z > = M%w*j> = -\dM\Hl\ 2 + |^°| 2 ] , (A.31) 

Mtw^z' = M B ] W3>Z , = gigz>[QH u \H° u \ 2 - Q Hd \H° d \ 2 } , (A.32) 

M B 2 m , z ,= 1 -g!i\H°\ 2 + \H°\ 2 }, (A.33) 

= M B 2 iW3>zl = g 2 gz'[QH d \H° d \ 2 - Q Hu \H° u \ 2 ] , (A.34) 

Mfw =2^(Q||^| 2 + g|J^| 2 + Q|jF d T + E^I^I 2 ) • (A.35) 

1=1 

In addition, for a scalar 0, we define 

A = [T 3 * - Q EM sin 2 %]M§ cos(2/3) , (A.36) 

3 

A; = ^[OjJflSl 2 + Q Hd \H° d \ 2 + Q S \S\ 2 + J2 QihW] ■ (A.37) 

i=l 

The up-type squark mass matrix is 

/ M 2 i h* u (hSH° d - At H°*) \ 

~ { KAh*S*H°* - A Ki Hl) Ml R ) ' {A - 6S) 

where 

M k = m k + < + A ^ + A k > (A.39) 

M% r = m% R + ml + A, k + A^ , (A.40) 

while the down-type matrix is obtained by u «-> d. The mass matrices for the sleptons 
and exotic scalars are similar. 
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B 11$ (T) for Scalar Particles in Model I 

The II$(T)s for the Higgs doublets are 

U Hd (T) = ^T 2 + l -giT> + \g 2 z ,Q\T 2 + £ (h% + \h 2 \ T 2 + \h 2 T 2 , (B.l) 

i=i ^ ' 



IWT) = \g 2 2 T 2 + i^ 2 T 2 + \g\,Q\T 2 + £ ^<T 2 + \h 2 T 2 . (B.2) 



i=i 

For the SM singlets: 

3 



n 5 (T) = ^I^r 2 + ia 2 T 2 + ^ 2 r 2 + ^^(af) 2 r 2 

i=l 

+ E^(«f) 2 ^ 2 + E7(«f)^ 2 > (B-3) 



2 v « , 
fc=l i=l 



n 5l (T) = ^|,Q|T 2 + iA 2 T 2 . (B.4) 

For the longitudinal components of the SU(3) C x SU{2) L x E/(l)y x C/(l)' gauge 
bosons, (7, W, B and Z': 

n 9 ,(T) = 6g 2 3 T 2 , U Wl (T) = 6g 2 2 T 2 , (B.5) 
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U Bl (T) = 8glT 2 , (T) = — g 2 z ,T 2 . (B.6) 

For Ni and S^: 

% i (T) = ^Q 2 v T 2 + i(af) 2 T 2 , (B.7) 



^(T)= 1 -g 2 z ,Q 2 Si T 2 + 1 -^ 



For X and X* 



n*(T) = ^i,g^r 2 + i« 2 T 2 , (b.9) 
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HxaCn = ^9 2 z'Qx 3 T 2 + -o?T 2 . (B.10) 

For Di and Df. 

n^(T) = \g\T 2 + l^T 2 + \&,<&F + \{afYT> , (B.ll) 

IL bi (T) = \g\T 2 + l^T 2 + ±^7* + \{affT 2 . (B.12) 
For i?; fc and H' d k : 

= ^ 2 2 T 2 + l^T 2 + ^ig^T 2 + ^(af ) 2 T 2 , (B.13) 

n^(T) = ^ 2 T 2 + i, 2 T 2 + ^I,g^T 2 + i(«f ) 2 T 2 . (B.14) 
For ^ and ^: 

n^(T) = ^ 2 T 2 + ^ 2 T 2 + ^g^T 2 , (B.15) 

n^(T) = ^ 2 T 2 + ^ 2 T 2 + ^I,g|,T 2 . (B.16) 
For the superpartners of the SM fermions: 

= \g\T 2 + \g\T* + 1, 2 T 2 + ^I,g Ql 2 T 2 + 1 « + /* 2 J T 2 , (B.17) 

n, k (T) = ^ 2 T 2 + 2 -g\T 2 + l&QvfT* + i<T 2 , (B.18) 

n 4 (T) = ^ 2 T 2 + ^gfT 2 + \&,Q3?T* + \h\T 2 , (B.19) 

= \g\T 2 + ^?T 2 + ^ 2 ,g Ls 2 T 2 + \h%T 2 , (B.20) 

n e - k (T) = \g\T 2 + i^g^T 8 + ^ 2 T 2 . (B.21) 
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C Derivation of AR in the Thick Wall Regime 



AR is the reflectivity asymmetry between the particle and antiparticle scatterings. It 
is crucial for the calculation of the CP- violating current in the EWBG. In the follow- 
ing, we will show how AR is calculated in the thick wall regime. We will explicitly 
calculate Ai? l y c ^^ c (/c i , kf, z) (see Eq. (BSD ) as an illustration. The calculations for 
ARt R _^t L (ki, kf, z) and AR tR _ >tL (ki, kf, z) (see Eq. flH7|) and Eq. fl99|) ) are similar, so 
we will not present their details. 

The general formula for the probability calculation in quantum field theory is 
(eg. see |131j ) 

d 3 p f 1 f d 3 pi ipi(pi) f d 3 pi ip*(pi) , . 



(2n) 3 2E f J {2nf^2Ej (2*)*^ 



Here the subscripts "i" and "f" denote initial and final states, respectively. Given 
that in the model the Langrangian is approximately invariant under the translations 
in the x, y and t directions, it is natural to define 



(P/bi) = (P/l* J d A xC(z)\pi) 

= t M(p l ^ Pf )(27r) 3 5(E t -E f )5(p*-p*)5(pV-p y f ), (C.2) 

with 



M = (pj| | d*£(*)|pf). (C3) 
Embedding Eq. IjCTl into Eq. flCHI we find 

where fc, i? and w are particle momentum, energy and velocity, respectively. For a 
given particle with mass m, all of the quantities are functions of k z . 

For the scattering process W c — > H c in the chargino sector, M(kf, kf, z) is 

M(k z , k), z) = j + ^ } dzi exp(^ 2 - kf)zi) 

(A(k z ,k z f ,z 1 )-A(k z ,k z f ,z 1 ^z)), (C.5) 

with 

A(k z ,k Z f,z x ) = UHo{k z f )[{u 1 {z 1 )e^{-id 1 {z l )-u l {z)ex^{-ie 1 {z))P L 

+(u 2 (zi) exp(i0 2 (zi) - u 2 (z) exp(i9 2 (z))P R ]u^ c (k z ), (C.6) 
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where Pl,r denote chirality projection operators. Then, 

A\M(k z ,k z f ,z)\ 2 = \M(k z ,k z f , z)\ 2 - (CP Conjugate) 



AA(k z ,k} , Z\ — > Z, Z 2 — > Z) 

-AA(k z , k z f , zi -> z, z 2 ) - AA(k?, k z f , z u z 2 ^ z), (C.7) 



with 



dz x dz 2 exp(i(kf - k z f )(z x - z 2 )) 
[A(k z , k z f , Zl )A(k z , k z f , z 2 )* - (CP Conjugate)]. (C.8) 

For the other AAs, only u Zi" or "z 2 " or both in the second line of Eq. (1C.8P are 
replaced by "2" . It is easy to calculate 

A.A(Af, Zi, z 2 ) = / dzxdz 2 sin((k z - k z f )( Zl - z 2 )) 

{fef^[ui(^i)ui(a? 2 ) sin(0i(*i) - 0i(z 2 )) 
+u 2 (zi)u 2 (z 2 ) sin(-0 2 (zi) + e 2 (z 2 ))\ + 
M 2; u(z)[ui(zi)w 2 (z 2 ) sin(0i(zi) + 9 2 (z 2 ) + a(z)) 
W^Kfe) sin(-0 2 (z 1 ) - 9 1 (z 2 ) - a(z))}}, (C.9) 

which implies 

AA(k z ,k z f , Zl ^z,z 2 ->z) = 0. (CIO) 

Finally, embedding Eq. (1C.9I) into Eq. (1C.7I) . we obtain Eq. (|93j) according to the 
relation 

AR^ c ^^ c (ki, k f , z) = -^-±-A\M(k z ,k z f ,z)\ 2 . (C.ll) 

4 A^/ v i v } 



D More Results at T = K in Model I 

Additional T = OA" results in the sMSSM with both ECPV and SCPV turned on 
are presented in Figure [16] - Figure [201 These results include the 7 dependences of 
tan p, m x c and m x o in cases a, b and c, and the correlations between tan (3, m x c and 
m o due to the variance of m 2 . For discussions on the phenomenology in a similar 
framework, also see [132] . 
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Figure 16: tan j3 vs. 7(71"). 
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Figure 17: m x c vs. 7(71"). 
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Figure 20: m x c vs. Tn x o. 
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